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Abstract 

It is shown in detail that the dynamics of the Einstein-dilaton-p- 
form system in the vicinity of a spacehke singularity can be asymp- 
totically described, at a generic spatial point, as a billiard motion in a 
region of Lobachevskii space (realized as an hyperboloid in the space 
of logarithmic scale factors) . This is done within the Hamiltonian for- 
malism, and for an arbitrary number of spacetime dimensions D > A. 
A key role in the derivation is played by the Iwasawa decomposition 
of the spatial metric, and by the fact that the off-diagonal degrees 
of freedom, as well as the p-form degrees of freedom, get "asymp- 
totically frozen" in this description. For those models admitting a 
Kac-Moody theoretic interpretation of the billiard dynamics we out- 
line how to set up an asymptotically equivalent description in terms 
of a one-dimensional non-linear cr-model formally invariant under the 
corresponding Kac-Moody group. 



1 Introduction 



1.1 BKL analysis in spacetime dimension D = 4 

The non-linearities of the Einstein equations are notably known to prevent 
the construction of an exact general solution. Only peculiar solutions, corre- 
sponding to idealized situations, have been explicitly derived. The singularity 
theorems of [48] predict the generic appearance of spacetime singularities un- 
der certain conditions, but do not provide a detailed description of how the 
spacetime becomes singular. From this perspective, the work of Belinskii, 
Khalatnikov and Lifshitz [6, 8], also known as "BKL", is quite remarkable as 
it gives a description of the generic asymptotic behaviour of the gravitational 
field in four spacetime dimensions in the vicinity of a spacelike singularity. 
As argued by these authors, near the singularity the spatial points essentially 
decouple, in the sense that the dynamical evolution of the spatial metric at 
each spatial point is asymptotically governed by a set of second-order, non- 
linear ordinary differential equations with respect to time. 

These differential equations are the same, at each spatial point, as those 
that arise in some spatially homogeneous cosmological models, which there- 
fore provide valuable insight into the qualitative features of the general so- 
lution. In the vacuum case, the spatially homogeneous models that capture 
the behaviour of the general solution are of Bianchi type IX or VIII (with 
homogeneity groups SU{2) or 5'L(2,M)). The asymptotic evolution of the 
metric can then be pictured as an infinite sequence of "oscillations" of the 
scale factors along independent spatial directions [6, 8]. This regime is called 
"oscillatory", or "of mixmaster type" [82], or also "chaotic" because it ex- 
hibits strong chaotic features [77, 16]. The coupling to matter fields does not 
change the picture, except if one includes a massless scalar field (equivalent 
to a perfect fiuid with "stiff' equation of state p — p), in which case the 
chaotic evolution is replaced by a monotonic power law evolution of the scale 
factors [5, 2], which mimicks the Kasner solution at each spatial point and is 
therefore called "Kasner-like" . In this case the spatially homogeneous model 
that captures the behaviour of the general solution is the Bianchi I model 
(with the abelian group of translations in as homogeneity group) . 
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1.2 BKL analysis in spacetime dimensions D > 4 



The extension of the BKL analysis to higher dimensions was addressed within 
the context of pure gravity (with no symmetry assumption) in [33, 32] , where 
it was shown that the general BKL approach remains valid: spatial points 
decouple as one approaches a spacelike singularity, i.e., the dynamical evolu- 
tion at each spatial point of the scale factors is again governed by ordinary 
differential equations. The main result of [33, 32], was that, while the general 
behaviour of solutions of the vacuum Einstein equations remains oscillatory 
for spacetime dimensions D < 10, it ceases to be so for spacetime dimensions 
D > 11, where it becomes Kasner-like. Let us also note that, just as in four 
spacetime dimensions, the coupling to a massless scalar field suppresses the 
chaotic behaviour in any number of spacetime dimensions and makes the 
solution monotonic (see, e.g., [29]). 

The authors of [33] did not consider the inclusion of massless p-iorms, 
which are part of the low energy bosonic sector of superstring/M-theory mod- 
els. This task was undertaken in [24, 25], with the finding that these p-forms 
play a crucial role and can reinstate chaos when it is otherwise suppressed. 
In particular, even though pure gravity is non-chaotic in eleven spacetime 
dimensions, the 3- form of D — 11 supergravity renders the system chaotic. 
Similarly, the bosonic sectors of all D = 10 supergravities related to string 
models define chaotic dynamical systems, thanks again to the p-forms, and 
in spite of the presence of a massless scalar dilaton. It is remarkable and sig- 
nificant that the (maximally supersymmetric) candidate models for a unified 
description of the fundamental forces not only have difficulties accomodating 
de Sitter-type spacetimes in any straightforward fashion [45, 80], but fur- 
thermore, and without exception, exhibit BKL chaos as one approaches the 
initial singularity. 

1.3 Billiard description of BKL behaviour 

An efficient way to grasp the asymptotic behaviour of the fields as one ap- 
proaches a spacelike singularity is based on the qualitative Hamiltonian meth- 
ods initiated by Misner [83] in the context of the Bianchi IX models (in four 
spacetime dimensions). The Hamiltonian approach naturally leads to a bil- 
liard description of the asymptotic evolution, in which the logarithms of the 
spatial scale factors define (after projecting out the dynamics of the overall 
volume factor) a geodesic motion in a region of the Lobachevskii plane H^, 
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interrupted by geometric reflections against the walls bounding this region 
[17, 84]. Chaos follows from the fact that the Bianchi IX biUiard has finite 
volume^. 

As pointed out in [71, 73, 56, 57, 26] this useful billiard description is quite 
general and can be extended to higher spacetime dimensions, with p-forms 
and dilaton. If d = D — 1 is the number of spatial dimensions, and if there are 
n dilatons, the billiard is a region of hyperbolic space Ha+n-i, each dilaton 
being equivalent, in the Hamiltonian, to the logarithm of a new scale factor. 
Besides the dimension of the hyperbolic billiard, the other ingredients that 
enter its definition are the walls that bound it. These walls can be of different 
types [24, 26] : symmetry walls related to the off-diagonal components of the 
spatial metric, gravitational walls related to the spatial curvature, and p- 
form walls (electric and magnetic) arising from the p-form energy- density. 
All these walls are hyperplanar. The billiard is a convex polyhedron with 
finitely many vertices, some of which are at infinity. In many important 
cases, the biUiard can be identified with the Weyl chamber of a Kac-Moody 
algebra, and the reflections against the billiard walls with the fundamental 
Weyl reflections [26, 27, 23]. This suggests deep connections with infinite 
symmetries. 

The main purpose of this paper is to provide a self-contained derivation 
of the billiard picture, in the general context of inhomogeneous solutions in 
D dimensions, with dilaton and p-form gauge fields. In particular, wc shall 
present a detailed derivation of the general results, announced and used in 
Ref. [24, 26], on the form of the various possible walls. For that purpose, 
we shall rely on the Iwasawa decomposition (see e.g. [49]) of the spatial 
metric. This provides an efficient derivation of the symmetry walls in any 
number of spacetime dimensions, which we obtain by working out explicitly 
the Hamiltonian that governs the dynamics in the "BKL limit" or "small 
volume limit". 

The present work is exclusively concerned with properties of solutions of 
Einstein's equations and their generalizations in the vicinity of a spacelike 
singularity. Accordingly, our treatment does not apply to timclike singular- 
ities, for which no analog of causal decoupling exists (the situation may be 
more subtle for the borderline case of a null singularity). Furthermore, we 

^Throughout this paper, the word biUiard used as a noun in the singular will denote 
the dynamical system consisting of a ball moving freely on a "table" (region in some 
Ricmannian space), with elastic bounces against the edges. Billiard will also sometimes 
mean the table itself. 
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are not making the claim here that all spacelike singularities are necessar- 
ily and uniformly of the BKL type. Rather, the results presented here are 
intended to refine and to generalize the work of [6, 7] by showing that, in 
higher dimensions and for many theoretically relevant matter sources, one 
can self- consistently describe the behaviour of all fields, in the vicinity of 
a spacelike singularity, and at a generic spatial point x, in terms of {%) a 
simple hyperbolic billiard description of the "angular" dynamics (7, vr^) of 
the logarithmic scale factors (after projection of the "radial" motion (p, vTp)), 
and {ii) an asymptotic "freezing" of the other phase space variables. This 
self-consistent asymptotic solution is general in the sense that it involves as 
many arbitrary functions of space as the most general solution. 

1.4 Organization of the paper 

After fixing our conventions and notations, and defining the class of La- 
grangians we shall consider in section 2, we discuss in section 3 the homo- 
geneous and diagonal Kasner solution in D spacetime dimensions, with a 
dilaton (as usual, the term "homogeneous" implies invariance under spa- 
tial translations). This solution plays a crucial role in the BKL approach 
because it describes the "free motion between collisions" and allows us to 
develop an important tool of our approach: the (Minkowskian) geometry of 
the scale factors. We then introduce (in section 4) the other main tool of 
our investigation: the Iwasawa decomposition of the spatial metric. To gain 
some familiarity with it, we discuss the asymptotic behaviour of non- diagonal 
Kasner metrics. 

In section 5, we explain the appearance of sharp potential walls in full 
generality, without imposing any homogeneity conditions on the metric and 
the matter fields. We then discuss in great detail the various walls that 
appear in physical models: symmetry walls, gravitational walls and p-form 
walls. The resulting geometry of the "billiard" made from all these walls 
is analyzed in section 7. Section 8 is devoted to the case when the biUiard 
can be identified with the Weyl chamber of a Kac-Moody algebra. To deal 
with this case we set up a Kac-Moody theoretic formulation of the billiard in 
terms of a non-linear cr-model based on the relevant Kac-Moody group in the 
last section. As we will show there, the asymptotic limit of these cr-models 
coincides with the asymptotic limit of the models discussed in the main 
part of this paper. An appendix illustrating by a toy model the asymptotic 
freezing of the off-diagonal degrees of freedom concludes this paper. 
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We should stress that our analysis is purely classical and accordingly, as 
it stands, is valid only up to the Planck (or string) scale. We shall also ignore 
fermionic fields throughout. Nevertheless, it is reasonable to expect that some 
of the ideas discussed here will remain relevant in a more general quantum 
mechanical context, at least qualitatively (see e.g. [88] for some recent ideas 
in this direction). The subject of Hamiltonian cosmology has a long history 
in the context of four-dimensional, spatially homogeneous spacetimes and 
provides useful insight on the general discussion presented here. For reviews 
on this subject, with an extensive bibliography, see [93, 94, 61]; see also the 
topical review on multidimensional gravity [59]. 

2 Models and Conventions 

2.1 The models 

We consider models of the general form 



R - dM(t>d^(t> 



4e 



2^(p + l)!" "A^i-Mp+i 



+ ... (2.1) 



where units are chosen such that IGttG/v = 1 (where is Newton's con- 
stant) and the spacetime dimension D = d + 1 is left unspecified. Besides 
the standard Einstcin-Hilbert term the above Lagrangian contains a dilaton 
field (p and a number of ]9-form fields A^^,,,j^^ (for p > 0). As a convenient 
common formulation we adopt the Einstein conformal frame and normalize 
the kinetic term of the dilaton with weight one w.r.t. to the Ricci scalar. 
The Einstein metric Gmn has Lorentz signature ( — !-■■■+) and is used to 
lower or raise the indices; its determinant is denoted by G. The p-form field 
strengths F^f) = dA^P^ arc normalized as 

^£J...M,+i = (P+l)^[Mi^Mt-M,+i] = ^MiAS^t-M,+i±P permutations . (2.2) 

The dots in the action (2.1) indicate possible modifications of the field 
strength by additional Yang- Mills or Chapline-Manton-type couplings [11, 
15], such as Fc = dG^'^^ - G^^^dB'^^^ for two 2-forms C^^) and E^^) and a 
0-form C*^°\ as they occur in type JIB supergravity. Further modifications 
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include Chern-Simons terms, as in the action for D = 11 supergravity [22]. 
The real parameter Ap measures the strength of the coupling of A'^^^ to the 
dilaton. When p = 0, we assume that Aq 7^ so that there is only one dilaton. 
This is done mostly for notational convenience. If there were other dilatons 
among the 0-forms, these should be separated off from the p-forms because 
they play a distinct role. They would define additional spacelike directions 
in the space of the (logarithmic) scale factors and would correspondingly 
increase the dimension of the relevant hyperbolic billiard. 

The metric Gmn, the dilaton field(s) and the p-form fields s-^® 
a priori arbitrary functions of both space and time, on which no symme- 
try conditions are imposed. Nevertheless it will turn out that the evolution 
equations near the singularity will be asymptotically the same as those of cer- 
tain homogeneous cosmological models. It is important to keep in mind that 
this simplification does not follow from imposing extra dimensional reduction 
conditions but emerges as a direct consequence of the general dynamics. 

2.2 Gauge conditions 

Our analysis applies both to past and future singularities, and in particular 
to Schwarzschild-type singularities inside black holes. To follow historical 
usage, we shall assume for definiteness that the spacelike singularity lies in 
the past, at finite distance in proper time. More specifically, we shall adopt 
a space-time slicing such that the singularity "occurs" on a constant time 
slice (t = in proper time). The slicing is built by use of pseudo- Gaussian 
coordinates defined by vanishing lapse N'^ — 0, with metric 

ds^ = -{N{x\ x')dx^f + gij{x^, x')dx'dx^ (2.3) 

In order to simplify various formulas later, we shall find it useful to introduce 
a rescaled lapse function 

N = N/^ (2.4) 

where g = deigij. We shall see that a useful gauge, within the Hamiltonian 
approach, is that defined by requiring 

N = p^ (2.5) 

where is a quadratic combination of the logarithms of the scale factors and 
the dilaton(s), which we will define below in terms of the Iwasawa decompo- 
sition. After fixing the time zero hypersurface the only coordinate freedom 
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left in the pseudo-Gaussian gauge (2.5) is that of making time-independent 
changes of spatial coordinates ,x* x'* = f \x^). Since the gauge condi- 
tion Eq.(2.5) is not invariant under spatial coordinate transformations, such 
changes of coordinates have the unusual feature of also changing the slicing. 
Throughout this paper, we will reserve the label t for the proper time 

dt = -Ndx^ = -N^dx\ (2.6) 

whereas the time coordinate associated with the special gauge Eq.(2.5) will 
be designated by T, viz. 

dT^-^. (2.7) 

Sometimes, it will also be useful to introduce the "intermediate" time coor- 
dinate r that would correspond to the gauge condition N — 1. It is explicitly 
defined by: 

dT^-—^p''dT (2.8) 

At the singularity the proper time t is assumed to remain finite and to de- 
crease toward 0+. By contrast, the coordinates T and r both increase toward 
-|-oo, as ensured by the minus sign in (2.6). Irrespective of the choice of co- 
ordinates, the spatial volume density g is assumed to collapse to zero at each 
spatial point in this limit. 

As for the p-form fields, we shall assume, throughout this paper, a gen- 
eralized temporal gauge, viz. 



= (2.9) 



where small Latin letter denote spatial indices from now on. This 

choice leaves the freedom of performing time-independent gauge transforma- 
tions, and therefore fixes the gauge only partially. 



3 Geometry of the space of the scale factors 

3.1 Supermetric and Hamiltonian 

To set the stage for our general Hamiltonian approach, it is useful to study 
first in detail the dynamics defined by considering only the kinetic terms 
of the metric and of the dilaton(s). The corresponding reduced action is 
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obtained from (2.1) by setting A^^ — and assuming that all the other fields 
depend only on time. In terms of a general time coordinate this reduced 
action reads 



S[g,„<P,N]= / dx^'N 



J(tr(g-ig)2-(trg-ig)2)+02 



(3.10) 



where we have suppressed an integral J d^x over the spatial volume for no- 
tational simplicity. Furthermore, we make use of the notations introduced 
in (2.3) and (2.4) with F = dF/dx^, and adopt a matrix notation where 
g{t) G GL{d,M^) stands for the matrix (gij) representing the spatial compo- 
nents of the metric at each spatial point. 

The action (3.10) is the (quadratic-in-velocities) action for a massless free 
particle with coordinates {gij, (p) moving in a curved target space with metric 

da' = \ [tr {g-Hgf - {irg-Hgf] + #2 (3.11) 

We designate by da'^ the line clement in this target "superspace" to distin- 
guish it from the line element in physical space time, which we denote by 
ds'. The first two terms in the r.h.s. of (3.11) define the so-called DeWitt 
supermetric in the space of the metric coefficients gij [34] . If several dilatons 
(ff (for i — 1, . . . ,n) were present the term dcf)^ in Eq.(3.11) would be re- 
placed by T,i{d(f)'^y. That is, each dilaton adds a (flat) direction in the target 
superspace. 

The rescaled lapse plays the role of an "einbein" in the geodesic action 
(3.10). As usual, extremizing over N yields the "zero-mass constraint" 

|(tr(g-^g)^-(trg-ig)2)+02^O. (3.12) 

Thus, the motion is given by a null geodesic of the metric (3.11). An affine 
parameter along those geodesies is dr — +Ndx° — —dt/y/g, cf. Eq.(2.8) 
above. In terms of the parameter r the equations of motion read: 

For diagonal metrics 

gK = exp [diag(-2/3)] ^ g^ = eM-^^ij (3-14) 
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the supermetric (3.11) reduces to 



d d 



Y^W? - (J] 0f + d(l>'' = G^,d(5''d(5r (3.15) 



1=1 i=i 



Here we have introduced a (d + 1) -dimensional space with coordinates (5'^^ 
with indices running over n = 1, ■ ■ ■ ,d + 1, such that the first d coordinates 
correspond to the logarithms of the scale factors of the spatial metric [cf. 
Eq. (3.14)], and the (d + l)-th coordinate j3'^^^ = represents the dilaton^. 
The explicit form of the (flat) target space metric Gjj^^ can be read ofT directly 
from (3.15). The action for diagonal metrics is 



In the sequel, we shall refer to this space as the "extended space of (logarith- 
mic) scale factors" or just "the /3-space" for short. 

Combining the scale factors and the dilaton(s) in a single space is natural 
because we know from Kaluza-Klein theory that the dilaton can be viewed as 
the logarithm of a scale factor in one extra spatial dimension. Independently 
of whether the original metric (2.3) has non-vanishing curvature or not, the 
metric (3.15) induced in the space of the scale factors (possibly including the 
dilaton) is flat. More precisely, the metric G^^, of the (extended) space of scale 
factors is a Minkowski metric in W^'^^ with signature ( — h H h). For in- 
stance, the direction in which only the dflaton varies [i.e. dfS^ oc (0, ■ ■ ■ , 0, 1)] 
is spacelike, while the direction in which only one scale factor varies [e.g., 
dp^ oc (1, 0, • • • , 0)] is null. A timelike direction in this space is the direction 
dl3^^ oc (1, 1, • • • , 1, 0). This reflects the famihar fact that the gravitational ac- 
tion is not bounded from below (even with Euclidean signature): conformal 
transformations of the metric, in which the scale factors arc all scaled in the 
same fashion, make da"^ negative. It is this characteristic feature of gravity 
which is responsible for the Lorentzian nature of the Kac-Moody algebras 
which emerge in the analysis of the billiard symmetries [26]. The Lorentzian 
signature of the metric in the space of the scale factors enables one to deflne 
the light cone through any point. We deflne the time-orientation to be such 

^Obviously, the range of indices would be extended io fi,y = 1, • • • , d+n in the presence 
of n dilatons. 




(3.16) 
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that future-pointing vectors v^^ have > 0. Geometrically, small volumes 

(small g) are associated with large positive values of Xlj/^*- Large volumes 
(large g), on the other hand, mean large negative values of We are 

interested in the small volume limit, i.e. —>■ +00. 

The Hamiltonian form of the action for the diagonal metric degrees of 
freedom and the dilaton is 

SIP'', TT^, N] = J [tt^/?'^ - INC^'tt^ttJ^ (3.17) 
where is the inverse of Gui,. Explicitly 

G'^%7r.^X^7rf-^('X:vr.) +4 (3.18) 

1=1 \i=l / 

where tt^ = (tTj, tt^) are the momenta conjugate to /3* and 0, respectively, i.e. 

TT^ = 2Ar-iG'^,/3- = 2(7^,-^ = 2(7^,^;'^ (3.19) 

or 

Here, the r-parameter velocities have been designated by f ^ = d/?^/ dr. 
3.2 Diagonal Kasner solution 

The Kasner solution (with or without dilaton) is now easily obtained by 
solving (3.12) and (3.13) in the diagonal case. Indeed, the equations of 
motion reduce to 

^ = 0. (3.20) 

They are solved by 

= vI't + (3.21) 

where v'^ and /3q are constants of the motion. The "zero-mass constraint" 
becomes 

G^^v^v" = 0. (3.22) 

One can transform the simple affine parameter solution (3.21) into the usual 
Kasner solution expressed in terms of the proper time by integrating the 
relation dt — —-^dr, with ^/g — exp(—Ej/?*), whence t oc exp ( — (EjV*)T), 
or ^ 

T — —= — r In t + const. (3.23) 
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We need to require Ej^;* > to remain consistent with our convention that 
T — > +00 at i — > 0"*" near the singularity. This yields 

d 

ds^ = -df + "^A^i{t){dx'f , Ai{t)^bitP' (3.24) 

1=1 

= -p^lnt + C^ (3.25) 

where bi = exp(— /3o) and = are integration constants and the minus 
sign in front of in (3.25) is included for the sake of uniformity in the 
formulas below (if there is no dilaton one simply sets = = 0). By 
rescaling the spatial coordinates, one can set bi — 1 and obtain the standard 
(proper time) form of the Kasner metric. The Kasner exponents = (Pi,P^) 
are given in terms of the affine velocities = dp^/ dr by 

P. = ^.- (3.26) 

Note that the sum in the denominator does not include the dilaton. 
They are subject to the quadratic constraint^ 

i=l \i=l ) 

coming from the "zero-mass condition" , and to the linear constraint 

Bi = 1 (3.28) 



coming from their definition above. 

If there are no dilatons, it follows from the above equations that there is 
at least one Kasner exponent which is negative, so at least one of the scale 
factors Aj(t) blows up as t ^ 0. The scale factors associated with positive 
Kasner exponents contract to zero monotonically. By contrast, in the pres- 
ence of a dilaton, all the Kasner exponents can be positive simultaneously. 
In both cases there is an overall contraction of the spatial volume since the 



■^Contrary to the variables and the velocities ^ we do not assign any covariance 
properties to the standard Kasner exponents, but regard them simply as parameters, 
leaving their labels always in the lower position. 
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determinant g of the spatial metric tends to zero. Indeed, a consequence of 
the hnear constraint above on the Kasner exponents is 

g oc t^. (3.29) 

Note that the relations (3.23) and (3.29) have been derived only for the exact 
(homogeneous) Kasner solution in the vacuum. 

3.3 Hyperbolic space 

Still as a preparation for dealing with the general inhomogeneous case, let 
us present an alternative way of solving the dynamics defined by the ki- 
netic terms of diagonal metrics, i.e. the action (3.16), or its Hamiltonian 
form. This alternative way will turn out to be very useful for describing 
the asymptotic dynamics of general inhomogeneous metrics. It consists in 
decomposing the motion of the variables {3^ into two pieces, namely a radial 
part p, and an angular one 7'^. Here "radial" and "angular" refer to polar 
coordinates in the Minkowski space of the More precisely, shifting, if 
necessary, the origin in /3-space to arrange that v^/3o < 0, the trajectories 
(3.21) will, for large enough values of r, get inside the future hght cone of 
the origin, i.e. 

= 2v^P^T + /?o ^/9o^ < 0. (3.30) 
Let us then decompose 13^ in hyperbolic polar coordinates (p, 7^*), i.e 

(3^ = P7^ (3.31) 

where 7^ are coordinates on the future sheet of the unit hyperboloid, which 
are constrained by 

I'l, = -1 (3.32) 

and p is the timelike variable 

p2 = -f3^^f3^ > (3.33) 

This decomposition naturally introduces the unit hyperboloid ("7-space"), 
which is a realization of the m-dimensional hyperbolic (Lobachevskii) space 
Hjn, with m — d — 1 + n if there are n > dilatons. 

In terms of the "polar" coordinates p and 7'', the metric in /3-space be- 
comes 

da^ = -dp"^ + p'^dE^ (3.34) 
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where dE^ is the metric on the 7-space if^. In these variables the Hamilto- 
nian reads 



Ho-- 



2 "^2 



(3.35) 



where tt^ are the momenta conjugate to the constrained hyperbolic coordi- 
nates 7'*. (It is straightforward to introduce angular coordinates on Hm to 
derive more explicit formulas for dY? and tt^.) The extra index on i7 is 
to underhne the fact that this Hamihonian refers only to a small part of the 
total Hamiltonian which we shall study below, because it describes only the 
kinetic terms associated to the {d + n) variables 13^. 
An equivalent expression is 

^o = ^[-7r^ + 7r^] (3.36) 

where we have introduced the new configuration variable 

A = lnp=|ln(-G'^,/3'^/3^) (3.37) 

with conjugate momentum tt^. 

The form (3.36) of the Hamiltonian shows that, when using a radial pro- 
jection on Hm-i the motion becomes simplest in the gauge 

N = (3.38) 

in terms of which (3.36) reduces to a free Hamiltonian on the pseudo-Riemann- 
ian space with metric —dX^ + In the gauge (3.38) the logarithmic radial 
momentum tta = pnp is a constant of the motion. In this gauge, we see that 
the free motion of the /3's is projected onto a geodesic motion on H^. 

The coordinate time T associated to this gauge, see Eq.(2.7), is hnked to 
the affine parameter r via 

dr 

dT^^. (3.39) 
P 

From (3.30) we get that varies linearly with r: 

p2 = -P^P'^ = -2v^P^T - P^Po^. (3.40) 

which implies 

r= ^lnT + const. (3.41) 
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Recalling that r varies logarithmically with the proper time t, we see that 
T oc In I lnt|. 

One can also use the configuration variable (3.37) as an intrinsic time 
variable to describe the dynamics. In view of (3.40) we have 

A = I In T + const. = | In | In i| + const. (3.42) 

Note that the various links above between the different time scales are derived 
only for an exact Kasner solution. Similar relations will hold asymptotically 
in the general inhomogeneous case (see discussion in chapter 7). 

4 Iwasawa decomposition and dynamics of 
non-diagonal Kasner metrics 

4.1 Iwasawa decomposition 

For homogeneous solutions in vacuo, the metric remains diagonal if the initial 
data are so. This is in general not true when matter (such as p-forms) or inho- 
mogeneities are included, in which case off-diagonal components generically 
appear even if there are none initially. For this reason, it is important to un- 
derstand the role of off-diagonal terms already in this simpler homogeneous 
context, by examining the evolution of initial data that are not diagonal. 
And, as for the simple diagonal metrics above, it is instructive to study the 
dynamics of non-diagonal metrics in several complementary ways. 

A first way of dealing with the evolution of initially non-diagonal metrics 
is to perform a suitable linear transformation on the diagonal Kasner solution. 
If L is the linear transformation needed to diagonalize the initial data, it is 
easy to see from (3.13) that the solution is given, in terms of the parameter 
T, by 

g(r) = L^g^(r)L (4.1) 

where ^ denotes transposition and gxiT) is the diagonal Kasner solution 
(3.14), (3.21). The dilaton, being a scalar, is still given by the same expression 
as before, i.e. = 13'^+^ = v'^+V + Note the relation 

detg(T) = (detL)2detg,^(T), (4.2) 

so that the relation (3.23) between r and t still holds. Therefore detg still 
goes to zero hke t'^. 



14 



A second way of describing the evolution of non-diagonal metrics is to 
perform an Iwasawa decomposition of the metric, 

g^M^A^N (4.3) 

where M is an upper triangular matrix with I's on the diagonal and ^ is a 
diagonal matrix with positive entries, which we parametrize as 

^ = exp(-/3), /3 = diag(/3^/3^•••,/3'^). (4.4) 

The explicit form of (4.3) reads 

d 
a=l 

The evolution of the new configuration variables = f{gij) defined by 
the Iwasawa decomposition differs in general from that of the quantities 
entering the diagonal Kasner solution (3.21), i.e. 7^ except for the 
special case jV^i = S^. Henceforth in this paper, the notation /3 will always 
refer to the logarithmic scale factors w.r.t. to the Iwasawa frame (4.5) (to 
be augmented, if needed, by the dilaton as a (rf + l)-th coordinate (3'^^^ = 
4>). Furthermore, we adopt the convention here and in the remainder that 
summation over the spatial coordinate indices is always understood, 

whereas sums over the Iwasawa frame indices a,b, . . . will always be written 
out. 

One can view the Iwasawa decomposition^ as the Gram-Schmidt orthog- 
onalisation of the initial coordinate coframe dx"^, which is indeed a triangular 
process, 

d 

gijdx'dx^ = e'' ® (4.6) 

a=l 

Starting with 9'^ = dx'^, one successively constructs the next 9's by adding lin- 
ear combinations of the dx'^ (for k > j) in such a way that O'''"^ is orthogonal 

*The Iwasawa decomposition applies to general symmetric spaces (see e.g. [49]). In 
our case the relevant symmetric space is the coset space SL{d,M.)/ SO{d) since the space 
of positive definite symmetric matrices can be identified with GL{d,'M.)/0{d), which is 
isomorphic to SL{d,R)/SOid) x R+. 
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to 9'^, 9'^ ^ is orthogonal to both 9*^ and 9'^ ^ etc. Explicitly, 



9' 



& 



id-l 



id-2 



= dx'^-^ +M'^-\dx'^, 



(4.7) 



Equivalently, 



= J\f% dx\ 



(4.8) 



Thus N'°'i vanishes for a > i, is equal to one for a = i, and has a non-trivial 
coordinate dependence only for a < i. It is easily seen that the determinant 
of the matrix J\f is equal to 1. Therefore the sum of the /5's is directly related 
to the metric determinant (as in the diagonal case): g — detg — det^^ = 



In the following we shall also need the vectorial frame {ea} dual to the 
coframe 9^: 



where the matrix jV*a is the inverse of A/'"j, i.e., M°'iN'^b — S^. It is again an 
upper triangular matrix with I's on the diagonal. 

4.2 Asymptotics of non-diagonal Kasner metrics 

Later in this paper we will deal with the generic inhomogeneous, non-diagonal 
metric gij{x^,x'^) and show how to study its Hamiltonian dynamics directly 
in terms of the Iwasawa variables {f5"'^M°'i) and of their conjugate momenta. 
But let us first take advantage of the known explicit solution (4.1) of the 
nondiagonal evolution in order to understand the qualitative behaviour of 
the Iwasawa variables (/?", AA^j) as t — > 0, i.e. r +oo. To be concrete let us 
explicitly consider the case of three spatial dimensions (c? = 3, D = 4) using 
the results of [51] for the Iwasawa decomposition of 3-dimensional metrics. 



exp(-2E„/?«). 




(4.9) 



Setting 



/ 1 ni n- 

1 
\0 1 




(4.10) 
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together with 



■exp(-/3i) 

\ exp(-/32) 

exp(-/?3^ 



one finds 



from which one gets 
n2 = 



-^In^ii, /3^ = -^ln 



9ll922 — 9l2 



-lln 

2 



9 



9ll922 — 9l2. 



Til 



9ii 

9l2 

9u 



9l3 923911 ~ 9l29l3 

— , ns = —. 

911 911922 - 9i2 



On the other hand, (4.1) and (3.24) yield 
with the constant matrix 

li I2 I3 

L = I mi 777.2 "^3 I , 



(4.11) 



9ii 




9i2 = riie 




9i3 ^n2e , 




922 






923 


= nin2e~'^'^ + n^e'^'^ , 




933 






+ e' 


-2/33 


(4.12) 



(4.13) 



(4.14) 



(4.15) 



ri r2 rz 

where wc have absorbed the ftj's in Zj, mi and r^. 

Combining these relations, we can deduce the explicit time dependence 
of the Iwasawa variables 



ni(t) 
n2{t) 



X 



-im 

2 



t"(^'+^'+p')(detL)- 
Y 



t^PHih + t^P^mim2 + t^P'rir2 
X ' 

X ' ^^^') = y 



(4.16) 
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with 

fP'+^P'{hm2 - l2mi){hm^ - hm,) 
^e^^+^V'^{hr2-l2ri){hn-hn) 

_l_^2p2+2p3 ^^^^2 — m2ri)(mir3 — msri). (4-17) 

These exphcit formulas show that the evolution of the Iwasawa variables 
(/?,jV) in the generic non-diagonal case is rather complicated. However, 
what will be important in the following is that they drastically simplify in 
the asymptotic hmit i ^ 0, i.e. r — > +00. Without loss of generahty, one 
can assume p\ < P2 < Pa- If necessary, this can be achieved by multiplying 
L by an appropriate permutation matrix. We shall in fact only consider 
Pi < P2 < Ps, leaving the discussion of the limiting cases to the reader. For 
generic L, i.e. h 7^ 0, ri 7^ 0, lim2 — Wmi 7^ 0, and mir2 — m2ri 7^ 0, one 
then finds the following asymptotic behaviour as i — > 0, i.e. r — > +00 : 



m = 

Y{t) = 
Z{t) = 



T — > +00 : ~ v^T, P"^ ~ V2T, 0^ ~ War, 

h h - Zgmi 

ni ^ n2^ -, na . (4.18) 

h li lim2 - hmi 

Here we used the links derived above between (p^, t) and {v^, r). 

The remarkable feature of these results is that, in the limit r +00, 
the evolution of the Iwasawa variables {P,J\f) is essentially as simple as the 
evolution we obtained in the diagonal case. Namely, the diagonal degrees of 
freedom /3 in the Iwasawa decomposition behave linearly with r, while the 
elements of the upper triangular matrix J\f tend to constants. That the n^'s 
should asymptotically tend to constants should be clear because they are 
homogeneous functions of degree zero in the metric coefficients — in fact, 
ratios of polynomials of degree one or two in the t^P\ It is more subtle that 
the scale factors exp(— 2/3*), which are homogeneous of degree one in the gi/s, 
are not all driven by the fastest growing (or least decreasing) term {t'^P^ for 
t — > 0+). This is what happens for the first scale factor exp(— 2/3^). However, 
the second scale factor exp(— 2/3^) feels the subleading term t^^'^ because the 
leading term drops from its numerator, equal to the minor giig2i — {gu)^- 
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Similar cancellations occur for the last scale factor exp(— 2/3^), which feels 
only the smallest term t^^^ as t — > 0^. 

The results (4.18) admit a simple generalization to d dimensions. By re- 
peating the above explicit calculation, one can prove that, in any dimension 
d, the /3's become asymptotically linear functions of r, as in the diagonal case, 
with coefficients that are given by a permutation of the underlying diagonal 
Kasner exponents. Furthermore, the A/'^j tend to constants, a phenomenon 
which we will refer to as the "asymptotic freezing" of the off-diagonal metric 
variables. More precisely, \ei < < ■ ■ ■ < v"^ be the ordered (unnormal- 
ized) underlying Kasner exponents; then, for generic L, we have 

/^"-'uV and AA^^ const. (4.19) 

as r — >^ +00 or, equivalently, t O"*". An alternative derivation of these 
results will follow, as a particular case, from the general result we shall de- 
rive below concerning the qualitative evolution of generic, inhomogeneous, 
matter-driven solutions in Iwasawa variables (after projecting out the radial 
motion of the Iwasawa (5 variables) . 

Note also that, in the simple case of homogeneous nondiagonal metrics 
in vacuo (without curvature and/or matter) one can discuss not only the 
limit t — > 0, but also the limit t — > -|-oo (i.e r — oo). In this second 
limit, one finds again that the /3's become asymptotically linear functions of 
T, with coefficients that are given by a permutation of the underlying Kasner 
exponents, and that the M°'i tend to constants. More precisely, for any d, 

^ yi.d-a) J. ^ -^hiie the explicit results in o? 3 read 

T -oo : ~ P3T , ~ P2T , (3^ ~ piT, 

ni ^ — , 712 — , ris^ (4.20) 

ri ri mir2 — m2ri 

Note that in both hmits, one has f3^ < p"^ <■■■ < p'^. The fact that the /3's 
are ordered in this way is due to our choice of an upper triangular A/"; had we 
taken J\f to be lower triangular instead, these inequalities would have been 
reversed. 



5 Asymptotic dynamics in the general case 

Having warmed up with the simple diagonal and homogeneous Kasner so- 
lution in Section 3, and having introduced the Iwasawa decomposition of 
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the metric (Section 4), we are now ready to apply these techniques to the 
description of the asymptotic dynamics in the hmit t to the general in- 
homogeneous, matter-driven case. The main ingredients (already introduced 
above for the homogeneous solutions) of our study are: 

• use of the Hamiltonian formalism, 

• Iwasawa decomposition of the metric, i.e. gij — > (/3",A/'"i), 

• decomposition of 13^ = (/3", 4>) into radial (p) and angular (7^) parts, 
and 

• use of the pseudo-Gaussian gauge (2.5), i.e. of the time coordinate T 
as the evolution parameter. 

More explicitly, with the conventions already described before, we assume 
that in some spacetime patch, the metric is given by (2.3) (pseudo-Gaussian 
gauge), such that the local volume g collapses at each spatial point as —>■ 
+00, in such a way that the proper time t tends to O"*". We work in the 
Hamiltonian formalism, i.e. with first order evolution equations in the phase- 
space of the system. For instance, the gravitational degrees of freedom are 
initially described by the metric Qij and its conjugate momentum vr*-'. We 
systematically use the Iwasawa decomposition (4.5) of the metric to replace 
the d{d + l)/2 variables g^j by the d + d{d - l)/2 variables (/3", A^"j). Note 
that {l3°',N"'i) are ultralocal functions of g^j, that is they depend, at each 
spacetime point, only on the value of gij at that point, not on its derivatives. 
This would not have been the case if we had used a "Kasner frame" (as 
defined below) instead of an Iwasawa one. The transformation g — > (/3,A/') 
then defines a corresponding transformation of the conjugate momenta, as we 
will explain below. We then augment the definition of the /3's by adding the 
dilaton field, i.e. = (/5",0), and define the hyperbolic radial coordinate p 
as in (3.33). Note that p is also an ultralocal function of the configuration 
variables {gij,(j)). We assume that the hyperbolic coordinate p can be used 
everywhere in a given region of space near the singularity as a well-defined 
(real) quantity which tends to -|-oo as we approach the singularity. We then 
define the slicing of spacetime by imposing the gauge condition (2.5). The 
time coordinate corresponding to this gauge is called T as above (see (2.7)). 
Our aim is to study the asymptotic behaviour of all the dynamical variables 
/3(T), A/"(T), as T — > -l-oo (recall that this limit also corresponds to i — > 0, 
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^ — > 0, p — > +00, with Z?'* going to infinity inside the future hght cone). 
Of course, we must also ascertain the self-consistency of this limit, which we 
shall refer to as the "BKL limit" . 



5.1 Hamiltonian action 



To focus on the features relevant to the billiard picture, we assume first that 
there are no Chern-Simons terms or couplings of the exterior form gauge 
fields through a modification of the curvatures F'^^\ which are thus taken 
to be Abelian, = dA^^ . We verify in subsection 6.5 below that these 
interaction terms do not change the analysis. The Hamiltonian action in any 
pseudo-Gaussian gauge, and in the temporal gauge (2.9), reads 



dx \ d X Ti'^g'ij + 7r00 + 



-y 



e \ (5.1) 



where the Hamiltonian density H is 
U = IC + M 



M 



1 i 1 1 2 



d-l 



~2p[ 



(5.2) 
(5.3) 



(5.4) 



-gR + gg^^d,cj>djcf> + ^ g F^.,^^^F^^^-=^-^ (5.5) 



where R is the spatial curvature scalar. The dynamical equations of motion 
are obtained by varying the above action w.r.t. the spatial metric compo- 
nents, the dilaton, the spatial p-form components and their conjugate mo- 
menta. In addition, there are constraints on the dynamical variables. 



n 

Hi 



("Hamiltonian constraint"), 
("momentum constraint"), 
( "Gauss law" for each p-form) 



(5.6) 
(5.7) 
(5.8) 
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with 



Hi 



p 
p 



= -^"-"u, (5-10) 

where the subscript \j stands for the spatially covariant derivative. 

Let us now see how the Hamiltonian action gets transformed when one 
performs, at each spatial point, the Iwasawa decomposition (4. 3), (4. 5) of 
the spatial metric. The "supermetric" (3.11) giving the kinetic terms of the 
metric and of the dilaton then becomes 

da^ = tTd(3' -{ti dp f + d(j)'' 

^ ^A\dJ\fJ\f-^)A-\dAfAf-Y] (5.11) 



I.e., 



da' = E(^^")' - [Y.'^f^j + 

a=l \a=l / 

_^^J2e'^f''-f'''^dAf%Af\Y (5.12) 

a<b 

where we recall that jV'a denotes, as in (4.9), the inverse of the triangular 
matrix J\f"'i appearing in the Iwasawa decomposition (4.5) of the spatial met- 
ric Qij. For d = 3, this expression reduces to the one of [51]. This change 
of variables corresponds to a point canonical transformation, which can be 
extended to the momenta in the standard way via 

a a 

Note that the momenta 



conjugate to the non-constant off-diagonal Iwasawa components J\f"-i are only 
defined for a < i; hence the second sum in (5.13) receives only contributions 
from a < i. 
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We next split the Hamiltonian density^ Ti (5.2) in two parts, one denoted 
by Tio, which is the kinetic term for the local scale factors (including 
dilatons) already encountered in section 3, and a "potential density" (of 
weight 2) denoted by V, which contains everything else. Our analysis below 
will show why it makes sense to group the kinetic terms of both the off- 
diagonal metric components and the p-forms with the usual potential terms, 
i.e. the term M. in (5.2). [Remembering that, in a Kaluza-Klcin reduction, 
the off-diagonal components of the metric in one dimension higher become 
a one-form, it is not surprising that it might be useful to group together the 
off-diagonal components and the p-forms.] Thus, we write 

n = no + V (5.15) 

with the kinetic term of the (5 variables 

Ho = ^G"^"7r^7r, (5.16) 

where the r.h.s. is that already defined in (3.18), with the replacement of 
the coordinate index i by the frame index a. The total (weight 2) potential 
density, 

V = V5 + Vg + ^Vp + V^, (5.17) 

p 

is naturally split into a centrifugal part linked to the kinetic energy of the off- 
diagonal components (the index "5"' referring to "symmetry", as discussed 
below) 

V5 = i5^e-^(^'-^")(P^•^^■)^ (5.18) 

a<h 

a "gravitational" (or "curvature") potential 

Vg = -gR. (5.19) 

and a term from the p-forms, 

V(P) = Vg) + V(7^" (5.20) 

^We use the term "Hamiltonian density" to denote both H and H. Note that is a 
usual spatial density (of weight 1, i.e. the same weight as y^), while H = ^/gH/N is a 
density of weight 2 (like g = (y^)^). Note also that tt*^ is of weight 1, while N = N/^ 
is of weight —1. 
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which is a sum of an "electric" and a "magnetic" contribution 



= (5-21) 



2p! 

ymagn _ jpip) ji-jp+i ( 99) 



Finally, there is a contribution to the potential linked to the spatial gradients 
of the dilaton: 

= gg''di(t>dj(t>. (5.23) 

We will analyze in detail these contributions to the potential, term by term, 
in section 6. 

5.2 Appearance of sharp walls in the BKL limit 
5.2.1 Derivation of central result 

In the decomposition of the Hamiltonian given above, we have split off the 
kinetic terms of the scale factors /3" and of the dilaton = from the 
other variables, and assigned the off-diagonal metric components and the p- 
form fields to various potentials, each of which is a complicated function of 
(3'^,N'"'i, P^a, ^jv--jp) ^(p) ''^ ^^'^ of some of their spatial gradients. The reason 
why this separation is useful is that, as we are going to show, in the BKL 
limit, and in the special Iwasawa decomposition which we have adopted, 
the asymptotic dynamics is governed by the scale factors whereas all 
other variables "freeze", just like for the general Kasner solution discussed 
in section 3. Thus, in the asymptotic limit, we have schematically 

V(r , AA^, P\, A^^),,^.^, 7^J;)■■^^ . . . ) ^ V^ir) (5.24) 

where Voo(7'^) stands for a sum of certain "sharp wall potentials" which 
depend only on the angular hyperbolic coordinates 7'' = p. As a con- 
sequence, the asymptotic dynamics can be described as a "billiard" in the 
hyperbolic space of the 7^'s, whose walls (or "cushions") are determined by 
the energy of the fields that are asymptotically frozen. 

This reduction of the complicated potential to a much simpler "effective 
potential" Vooil'^) follows essentially from the exponential dependence of 
V on the diagonal Iwasawa variables (3^, from its independence from the 
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conjugate momenta of the /5's, and from the fact that the radial magnitude 
p of the /3's becomes infinitely large in the BKL limit. 

To see the essence of this reduction, with a minimum of technical com- 
plications, let us consider a potential density (of weight 2) of the general 
form 

V(/3, g, P) = J] ca{Q, P) exp ( - 2wAm) (5-25) 

A 

where (Q, P) denote the remaining phase space variables (that is, other than 
(/9, VTa)). Here wa_{I3) = w^f^f]^ are certain linear forms which depend only 
on the (extended) scale factors, and whose precise form will be derived in 
the following section. Similarly we shall discuss below the explicit form of 
the pre-factors ca, which will be some complicated polynomial functions of 
the remaining fields, i.e. the off-diagonal components of the metric, the p- 
form fields and their respective conjugate momenta, and of some of their 
spatial gradients. The fact that the wa{P) depend linearly on the scale 
factors P'^ is an important property of the models under consideration. A 
second non-trivial fact is that, for the leading contributions, the pre-factors 
are always non-negative, i.e. q^'^^^^ > Q. Since the values of the fields for 
which Ca = constitute a set of measure zero, we will usually make the 
"genericity assumption" ca > for the leading terms in the potential V®. 
The third fact following from the detailed analysis of the walls that we shall 
exploit is that all the leading walls are timelike, i.e. their normal vectors (in 
the Minkowski /5-space) are spacelike. 

As shown in section 3.3, the part of the Hamiltonian describing the kinetic 
energy of the /5's, Hq = NTCq, takes the form (3.36) when parametrizing 
in terms of p and 7^, or equivalently, A = Inp and 7'' (cf. Eq. (3.37)). 
Choosing the gauge (2.5) to simplify the kinetic terms Hq we end up with an 
Hamiltonian of the form 

H{\,nx,^,n^,Q,P) = NH (5.26) 

A 

where tt^ is the kinetic energy of a particle moving on Hm- In (5.26) and 
below we shall regard A as a primary dynamical variable (so that p = e'^). 

^Understanding the eflFects of the possible failure of this assumption is one of the subtle 
issues in establishing a rigorous proof of the BKL picture. 
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The essential point is now that, in the BKL hmit, A +00 i.e. p — * +00, 
each term exp ( — 2pwA{l)) becomes a sharp wall potential, i.e. a function 
of wa{'^) which is zero when ^^(7) > 0, and +00 when wa{'^) < 0. To 
formalize this behaviour we define the sharp wall 0-function^ as 

A basic formal property of this ©-function is its invariance under multiplica- 
tion by a positive quantity. With the above assumption checked below that 
all the relevant prefactors ca{Q,P) are positive near each leading wall, we 
can formally write 



lim 

p— >oo 



ca{Q, P)p'exp ( - pwAi-f)] = ca{Q, P)©( - 2wA{-f)) = ©( - 2wA{-f)). 

(5.28) 

Of course, Q{—2wA{'y)) — ©(— wa(7)), but we shall keep the extra factor of 2 
to recall that the arguments of the exponentials, from which the ©-functions 
originate, come with a well-defined normalization. Therefore, the limiting 
Hamiltonian density reads 

HooiX, TTA, 7, ^y, Q,P) = \ [-^l + ^7] + E ®( - 2«'A'(7)), (5.29) 

A' 

where A' runs over the dominant walls. The set of dominant walls is defined 
as the minimal set of wall forms which suffice to define the billiard table, i.e. 
such that the restricted set of inequalities {wA'i'y) > 0} imply the full set 
{^^(7) > 0}. The concept of dominant wall will be illustrated below. [Note 
that the concept of "dominant" wall is a refinement of the distinction, which 
will also enter our discussion, between a leading wall and a subleading one.] 
The crucial point is that the limiting Hamiltonian (5.29) no longer de- 
pends on A, Q and P. Therefore the Hamiltonian equations of motion for 
A, Q and P tell us that the corresponding conjugate momenta, i.e. tt^, P and 
Q, respectively, all become constants of the motion in the limit A +00. 
The total Hamiltonian density is also a constant of the motion (which 
must be set to zero). The variable A evolves according to dX/dT — — ^tta. 
Hence, in the limit, A is a linear function of T. The only non-trivial dynamics 



^Onc should more properly write Ooo(2;), but since this is the only step function en- 
countered in this article, we use the simpler notation Q{x). 
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resides in the evolution of (7, n^) which reduces to the sum of a free (non 
relativistic) kinetic term and a sum of sharp wall potentials, such that 
the resulting motion of the 7's indeed constitutes a billiard, with geodesic 
motion on the unit hyperboloid Hm interrupted by reflections on the walls 
defined by WAij) = 0. These walls are hyperplanes (in the sense of hyper- 
bolic geometry) because they are geometrically given by the intersection of 
the unit hyperboloid (3^(3^ — —1 with the usual Minkowskian hyperplanes 



We note in passing that an alternative route for reducing the dynamics 
to a billiard in 7-space would be to eliminate the variable A by solving the 
Hamiltonian constraint for n^. This allows one to use A as a time variable, 
and is similar to going from the quadratic form of the action of a relativistic 
particle to its "square-root form" 



As the above derivation of the asymptotic constancy of all the phase-space 
variables {Q, P) may seem a bit formal, we study, in Appendix A, a simplified 
model explicitly showing how this asymptotic constancy arises. This toy 
model also exemplifies the residual, asymptotically decaying variations of 
(Q, P) as p — >• 00. 

5.2.2 Finite volume vs. infinite volume 

Geodesic motion in a billiard in hyperbolic space has been much studied. It 
is known that this motion is chaotic or non-chaotic depending on whether 
the billiard has finite or infinite volume [81, 52, 101, 36]. In the finite volume 
case, the generic evolution exhibits an infinite number of collisions with the 
walls with strong chaotic features ("oscillating behavior"). 

By contrast, if the billiard has infinite volume, the evolution is non- 
chaotic. For a generic evolution, there are only finitely many collisions with 
the walls. The system generically settles after a finite time in a Kasner-like 
motion that lasts all the way to the singularity. 



wa{P) = 0. 
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5.2.3 /3-space description 

The above derivation relied on the use of hyperbohc polar coordinates (p, 7). 
This use is technically useful in that it represents the walls as being located 
at an asymptotically fixed position in hyperbolic space, namely ^^(7) = 0. 
However, once one has derived the final result (5.29), one can reexpress it 
in terms of the original variables /S'^, which run over a linear (Minkowski) 
space. Owing to the linearity of wa{P) = pwAij) in this Minkowskian pic- 
ture, the asymptotic motion takes place in a "polywedge" , bounded by the 
hyperplanes wa{P) = 0. The billiard motion then consists of free motions 
of (3^ on straight lightlike lines within this polywedge, which are interrupted 
by specular reflections off the walls. [See formula (7.3) below for the explicit 
effect of these reflections on the components of the velocity vector of the 
/3-particle.] Indeed, when going back to /3-space (i.e. before taking the BKL 
limit), the dynamics of the scale factors at each point of space is given by 
the Hamiltonian 



iG'^'n^Ti, + J] CA exp ( - 2wa{P)) 



(5.31) 



The /3-space dynamics simplifies in the gauge N — 1, corresponding to the 
time coordinate r. In the BKL limit, the Hamiltonian (5.31) takes the lim- 
iting form (in the gauge N — 1) 

HUP", 71^) = iG'^%7r, + ^e(-2wA'(/3)) (5.32) 

A' 

where the sum is again only over the dominant walls. When taking equal 
time slices of this polywedge (e.g., slices on which is constant), it is 
clear that with increasing time (i.e increasing or increasing r or p) the 

walls recede from the observer. The /3-space picture is useful for simplifying 
the mathematical representation of the dynamics of the scale factors which 
takes place in a linear space. However, it is inconvenient both for proving 
that the exponential walls of (5.31) do reduce, in the large limit to 
sharp walls, and for dealing with the dynamics of the other phase-space 
variables {Q, P), whose appearance in the coefficients ca has been suppressed 
in (5.31) above. Let us only mention that, in order to prove, in this picture, 
the freezing of the phase space variables (Q, P) one must consider in detail 
the accumulation of the "redshifts" of the energy-momentum tt^ of the f3- 
particle when it undergoes reflections on the receding walls, and the effect 
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of the resulting decrease of the magnitudes of the components of tt^ on the 
evolution of {Q,P). (In the notation of Appendix A, it is important to take 
into account the fact that po decreases with time.) By contrast, the 7-space 
picture that we used above allows for a more streamlined treatment of the 
effect of the limit p — > +00 on the sharpening of the walls and on the freezing 

of (g,p). 

In summary, the dynamics simplifies enormously in the asymptotic limit. 
It becomes ultralocal in that it reduces to a continuous superposition of 
evolution systems (depending only on a time parameter) for the scale factors 
and the dilatons, at each spatial point, with asymptotic freezing of the off- 
diagonal and p-form variables. This ultralocal description of the dynamics 
is valid only asymptotically. It would make no sense to speak of a billiard 
motion prior to this limit, because one cannot replace the exponentials by 
0-functions. Prior to this limit, the evolution system for the scale factors 
involves the coefficients in front of the exponential terms, and the evolution 
of these coefficients depends on various spatial gradients of the other degrees 
of freedom. However, one may contemplate setting up an expansion in which 
the sharp wall model is replaced by a model with exponential ( "Toda-like" ) 
potentials, and where the evolution of the quantities entering the coefficients 
of these "Toda walls" is treated as a next to leading effect. See [28] for 
the definition of the first steps of such an expansion scheme for maximal 
supergravity in eleven dimensions. 

5.3 Kasner frames vs. Iwasawa frames 

At this point, we want to clarify an issue that might at first sight seem 
paradoxical to the reader. In most of the BKL literature (notably in the 
original analysis of [6, 8]), one tries to construct the metric in frames equal 
(or close) to "Kasner frames"^. These are defined as frames with respect 
to which both the spatial metric and the extrinsic curvature are diagonal. 
Once the time slicing has been fixed, this geometric frame is unique, up 
to time- dependent rescalings of each basis vector, when the eigenvalues of 
the extrinsic curvature are distinct and some definite ordering of the (time- 
independent) eigenvalues Pa has been adopted. It is not clear how to fix in 
a rigourous (and useful) manner the arbitrary time-dependent rescalings of 

^We know, however, no detailed development of the Hamiltonian formalism within such 
a frame whose definition involves both the metric variables and their conjugate momenta. 



29 



each basis vector. However, this can be done (following [7]) in an approximate 
manner by considering the evolution of the geometry at each spatial point 
as a succession of free Kasner flights interrupted by collisions against the 
walls. More precisely, one can uniquely fix the normalization of the Kasner 
frame by imposing two requirements: (i) that the Kasner frame be time- 
independent during each Kasner epoch, i.e. that the Lie derivative of the 
frame vectors along d/dt be zero, and (ii) that the linear transformation 
between two successive Kasner frames has the special form given in (5.35) 
below (and not (5.35) up to some time-independent rescaling). 

More explicitly, if we consider the three-dimensional case for definiteness 
and denote the covariant components of the frame by {Zi,mj,rj}, one has 
(suppressing the x-dependence) 

g,j{t) = Al{t)lilj + Al{t)m,mj + Al{t)r,rj (5.33) 

during a certain Kasner epoch, and 

g,,it) = Ajit) + Alit) rri^n-i'^ + Al{t) r'f'^ (5.34) 

during the subsequent Kasner epoch. Here the two successive Kasner frames 
{/i, mj, Tj}, {/^, m^, r[} are (in this approximation) independent of time, while 
the scale factors Aait) vary like a power law during each Kasner free flight (say 
Aa{t) ^ hatP" during the first epoch, and Aa{t) ^ h'Jp''^ during the next, with 
some interpolating behaviour during the collision). Ref. [7] has argued (by 
studying the effect of one collision in the "incoming" frame {/j,mj,rj}, and 
by rediagonalizing the "outgoing" metric) that the transformation between 
the two successive Kasner frames could be written as: 

= li, m- = mi GmK r\^ri + OrU (5.35) 

Here, cr^ and cr^ are quantities which can a priori be of order unity. 

The seeming paradox is that according to [7] the Om and Or do not get 
smaller for collisions closer to the singularity, so that the transformation 
(5.35) from the old Kasner axes / to the new Kasner axes /' is of order one, 
no matter how close one gets to the singularity; therefore, the Kasner axes 
generically never "come to rest" if there is an infinite number of collisions. On 
the other hand, as we just saw, the Iwasawa frames become approximately 
time-independent for asymptotic values of r, and the extrinsic curvature 
approximately diagonal. 
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We now show that there is no contradiction between the oscillatory behav- 
ior of the Kasner axes (in the generic inhomogeneous case) and the asymp- 
totic freezing of the off-diagonal components M in the Iwasawa frame. The 
key point is the very restricted form of the change (5.35) of the covariant 
components during a collision. Recall first that the covariant components of 
the Iwasawa frame are given the three covectors appearing as the three lines 
of the matrix A/", i.e. 9^ = dx^ + riidx^ + n2dx^, (P' = dx'^ + n^dx^, 6^ = dx^. 
[Note also in passing that, contrary to the Kasner frames, the Iwasawa frames 
are not geometrically uniquely defined since one can redefine the coordinates 
x\] 

Let us now see what (5.35) imphes for the change in the Iwasawa variables. 
For doing this we need to relate these to the rui, rii. This was done in 
section 4.2 for the simple exact Kasner solution and is easy to generalize to 
the case of an oscillatory metric. Indeed, when one is not in the "collision 
region" (and this can be applied both before and after the specific collision 
under consideration in (5.35)), the coordinate components of the metric take 
the form (5.33) or (5.34) with Al < A\ and Al < Al. As the definition of 
the Iwasawa components Ua is purely algebraic, we can apply the formulas 
(4.16)-(4.17) of section 4.2 to the present case. It suffices to replace t^" by 
Aa everywhere. Then, it is easy to see that ni, n2 and are still given by 
the same final formula as above, i.e 

ni = — , 712 = -r, ns = — j-— , (5.36) 

h li 

(before the collision) and 



^' - - 7i' - 7l' 





- hrrii 




— kmi ' 







(after the collision). If we substitute in this second formula m[ and r[ in 
terms of Zj, mi and n according to (5.35), we get the same values for the 
Iwasawa off-diagonal variables rii, n2 and ^3, before and after the collision, 
namely, n\ = ni, n'2 = n2, n'^ = n^. There is thus no contradiction between 
the change of Kasner axes (5.35) and the freezing of the off-diagonal Iwasawa 
variables. The same conclusion holds for collisions against the other types of 
walls, where the Kasner axes "rotate" as in (5.35). 
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5.4 Constraints 



We have just seen that in the BKL hmit, the evolution equations become 
ordinary differential equations with respect to time. Although the spatial 
points are decoupled in the evolution equations, they are, however, still cou- 
pled via the constraints. These constraints just restrict the initial data and 
need only be imposed at one time, since they are preserved by the dynamical 
equations of motion. Indeed, one easily finds that, in the BKL limit. 



since [H{x),H{x')] = in the ultralocal hmit. This corresponds simply to 
the fact that the collisions preserve the lightlike character of the velocity 
vector. Furthermore, the gauge constraints (5.8) are also preserved in time 
since the Hamiltonian constraint is gauge- invariant. In the BKL limit, the 
momentum constraint fulfills 



It is important that the restrictions on the initial data do not bring dangerous 
constraints on the coefficients of the walls in the sense that these may all take 
non-zero vahics. For instance, it is well known that it is consistent with the 
Gauss law to take non-vanishing electric and magnetic energy densities; thus 
the coefficients of the electric and magnetic walls are indeed generically non- 
vanishing. In fact, the constraints are essentially conditions on the spatial 
gradients of the variables entering the wall coefficients, not on these variables 
themselves. In some non-generic contexts, however, the constraints could 
force some of the wall coefficients to be zero; the corresponding walls would 
thus be absent. [E.g., for vacuum gravity in four dimensions, the momentum 
constraints for some Bianchi homogeneous models force some symmetry wall 
coefficients to vanish. But this is peculiar to the homogeneous case.] 

It is easy to see that the number of arbitrary physical functions involved 
in the solution of the asymptotic BKL equations of motion is the same as 
in the general solution of the complete Einstein-matter equations. Indeed, 
the number of constraints on the initial data and the residual gauge freedom 
are the same in both cases. Further discussion of the constraints in the BKL 
context may be found in [2, 29]. 



(5.38) 



(5.39) 
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5.5 Consistency of BKL behaviour in spite of the in- 
crease of spatial gradients 

The essential assumption in the BKL analysis, also made in the present pa- 
per, is the asymptotic dominance of time derivatives with respect to space 
derivatives near a spacelike singularity. This assumption has been mathe- 
matically justified, in a rigorous manner, in the cases where the billiard is of 
infinite volume, i.e. in the (simple) cases where the asymptotic behaviour is 
not chaotic, but is monotonically Kasner-like [2, 29]. 

On the other hand, one might a priori worry that this assumption is 
self-contradictory in those cases where the billiard is of finite volume, when 
the asymptotic behaviour is chaotic, with an infinite number of oscillations. 
Indeed, it has been pointed out [72, 4] that the independence of the billiard 
evolution at each spatial point will have the effect of infinitely increasing 
the spatial gradients of various quantities, notably of the local values of the 
Kasner exponents p^{x). This increase of spatial gradients towards the singu- 
larity has been described as a kind of turbulent behaviour of the gravitational 
field, in which energy is pumped into shorter and shorter length scales [72, 4], 
and, if it were too violent, it would certainly work against the validity of the 
BKL assumption of asymptotic dominance of time derivatives. For instance, 
in our analysis of gravitational walls in the following section, wc will en- 
counter subleading walls, whose prefactors depend on spatial gradients of 
the logarithmic scale factors /3. 

To address the question of consistency of the BKL assumption we need to 
know how fast the spatial gradients of /3, and of similar quantities entering 
the prefactors, grow near the singularity. Let us consider the spatial gradient 
oi (3 = fTf, which is 

dip = pda + PldiX (5.40) 

As we are working here in the gauge (2.5), the spatial derivatives must be 
taken with fixed T. We know that, asymptotically, A is a linear function of 
T, i.e. A = a{x)T + b{x) where a{x) — —^7rx{x) is linked to the (spatially 
dependent) conjugate momentum of A. Therefore the spatial gradient diX — 
Tdia{x) + dib{x) behaves linearly in T, so that diX oc A = In p. The second 
term in (5.40) consequently behaves as 7/7 In p when p 00. Let us now 
estimate the first term (which will turn out to dominate the sum). 

To estimate 9j7 we can use the standard results on biUiards on hyperbolic 
space. Indeed, dx^di'y can be thought of as the infinitesimal deviation between 
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two billiard trajectories. The chaotic behaviour of the billiard implies that 
this deviation will grow exponentially with T. In fact, because we work 
on Hm with curvature —1 the Liapunov exponent for billiard trajectories is 
equal to one, when using geodesic length as time parameter. Remembering 
the Hamiltonian constraint tt^ = tt^, the geodesic length is simply equal to 
A. This yields the growth estimate 

da = exp A = 0{l)p (5.41) 

where the coefficient 0{1) is a chaotically oscillating quantity. This estimate 
is not affected by the collisions against the walls because these preserve the 
angles made by neighboring trajectories (the walls are hyperplanes). Insert- 
ing (5.41) in (5.40), we see that the first term indeed dominates the second. 
We conclude that the chaotic character of the billiard indeed imphes an un- 
limited growth of the spatial gradients of (3, but that this growth is only of 
polynomial order in p 

d,(3 = 0{l)p^ (5.42) 

This polynomial growth of dif3 (and of its second-order spatial derivatives) 
entails a polynomial growth of the prefactors of the sub-dominant walls in 
section 6.2. Because it is polynomial (in p), this growth is, however, negligible 
compared to the exponential (in p) behaviour of the various potential terms. 
It does not jeopardize our reasoning based on keeping track of the various 
exponential behaviours. As we will see the potentially dangerously growing 
terms that we have controlled here appear only in subdominant walls. The 
reasoning of the Appendix shows that the prefactors of the dominant walls 
are self-consistently predicted to evolve very little near the singularity. 

We conclude that the unlimited growth of some of the spatial gradients 
does not affect the consistency of the BKL analysis done here. This does 
not mean, however, that it will be easy to promote our analysis to a rigorous 
mathematical proof. The main obstacle to such a proof appears to be the 
existence of exceptional points, where a prefactor of a dominant wall happens 
to vanish, or points where a subdominant wall happens to be comparable to 
a dominant one. Though the set of such exceptional points is (generically) 
of measure zero, their density might increase near the singularity because of 
the increasing spatial gradients. This situation might be compared to the 
KAM (Kolmogorov-Arnold-Moser) one, where the "bad" tori have a small 
measure, but are interspersed densely among the "good" ones. 
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5.6 BKL limit vs. other limits 



The last issue which we wish to address in this chapter concerns the relation 
of the BKL limit to other limits considered in the hterature, such as the 
strong coupling limit, or the "small tension limit", as well as the relation 
with asymptotically velocity dominated solutions. 

It is sometimes useful to separate the time derivatives (conjugate mo- 
menta) in the Hamiltonian from the space derivatives, viz. 

n^}C' + eV (5.43) 

where £ = ±1 according to whether the spacetime signature is Lorentzian 
{e — 1) or Euchdean {e — —1). Here, 

IC' = Ho + Vs + Vg) (5.44) 

contains all the kinetic terms, and 

V' = Vg + V™;^" + V^ (5.45) 

the terms with spatial derivatives. We stress that this split is different from 
the one introduced in subsection 5.1, where only the kinetic terms of the 
sclae factors and the dilatons were kept as such. The above split is useful 
because for some models the asymptotic dynamics is entirely controlled by 
/C', i.e., by the hmit £ = 0. This happens whenever the biUiard that emerges 
in the BKL limit is defined by the symmetry and electric walls, as for in- 
stance for D = 11 supergravity [24], or the pure Einstein-Maxwell system in 
spacetime dimensions D > 5 [25, 74]. Curvature and magnetic walls are then 
subdominant, i.e., spatial gradients become negligible as one approaches the 
singularity. 

If the curvature and magnetic walls can be neglected, the evolution equa- 
tions are exactly the same as the equations of motion obtained by performing 
a direct torus reduction to 1 -|- dimensions. We stress, however, that no ho- 
mogeneity assumption has to be made here. The effective torus dimensional 
reduction follows from the dynamics and is not imposed by hand. 

The limit £ = is known as the "zero signature limit" [97] and lies 
half-way between spacetimes of Minkowskian or Lorentzian signature. It 
corresponds to a vanishing velocity of light (or vanishing "medium tension"); 
the underlying geometry is built on the Carroll contraction of the Lorentz 
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group [50]. The terminology "strong coupling" is also used [55] and stems 
from the fact that, with appropriate rescalings, H. can be rewritten as 

H = GnK' + G-^'^V (5.46) 

such that the limit in question corresponds to large values of Newton's con- 
stant Gn- The interest in this ultrarelativistic ( "CarroUian" ) limit has re- 
cently been revived in [30, 78, 1] and [46]. 

If there are only finitely many collisions with the walls (corresponding to 
a billiard of infinite volume) the dynamics in the vicinity of the singularity 
becomes even simpler. After the last collision, the asymptotic dynamics is 
controlled solely by the kinetic energy TYq of the scale factors. This case where 
both spatial gradients and matter (here p-ioim) terms can be neglected, 
has been called "asymptotically velocity-dominated" (or "AVD" ) in [35] and 
allows a rigorous analysis of its asymptotic dynamics by means of Fuchsian 
techniques [2, 29]. By contrast, rigorous results are rare for the case of 
infinitely many collisions (sec, however, the recent analytic advances in [91]). 
Besides the existing rigorous results, there also exists a wealth of numerical 
support for the BKL ideas [10, 9]. 

6 Walls 

The decomposition (5.15) of the Hamiltonian gives rise to different types 
of walls, which we now discuss in turn. Specifically, we will derive explicit 
formulas for the linear forms wa{13) and the field dependence of the pre- 
factors ca entering the various potentials. 

6.1 Centrifugal (or symmetry) walls 

We start by analyzing the effects of the off-diagonal metric components which 
will give rise to the so-called "symmetry walls" . As they originate from the 
gravitational action they are always present. The relevant contributions to 
the potential is the centrifugal potential (5.18). When comparing (5.18) to 
the general form (5.25) analyzed above, we see that firstly the summation 
index A must be interpreted as a double index (a, b) , with the restriction 
a < b, secondly that the corresponding prefactor is c^b = {P^Af^jY is auto- 
matically non-negative (in accordance with our genericity assumptions, we 
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shall assume Cab > 0). The centrifugal wall forms read: 

^fa6)(/3) ^ wf,b),P'^ ^(5'-(5^ {a< h). (6.1) 

We refer to these wall forms as the "symmetry walls" for the following reason. 
When applying the general collision law (7.3) derived below to the case of the 
collision on the wall (6.1) one easily finds that its effect on the components of 
the velocity vector v'^ is simply to permute the components and v^, while 
leaving unchanged the other components /i ^ a,b. 

The hyperplanes w^ab-^iP) — (i.e. the symmetry walls) are timelike since 

G'"'y^iab),y^iab)u = +2 (6.2) 

This ensures that the symmetry walls intersect the hyperboloid Gn^P^P'^ — 
— 1, > 0. The symmetry billiard (in /3-space) is defined to be the 

region of Minkowski space determined by the inequalities 

«^fa6)(/3)>0, (6.3) 

with 0"" > (i.e. by the region of /3-space where the functions are 
zero). Its projection on the hyperbolic space is defined by the inequalities 

^fa6)(7) > 0. 

The explicit expressions above of the symmetry wall forms also allow us to 
illustrate the notion of a "dominant wall" defined in subsection 5.2.1 above. 
Indeed, the d{d — l)/2 inequalities (6.3) already follow from the following 
minimal set of ci — 1 inequalities 

- >0, - >0, - P'^'^ > (6.4) 

More precisely, each linear form which must be positive in (6.3) can be written 
as a linear combination, with positive (in fact, integer) coefficients of the 
linear forms entering the STibset (6.4). For instance, f3^ — f3^ = (/?'^ — /3^) + 
{j3'^ — j3^), etc. In the last section, we will reinterpret this result by identifying 
the dominant linear forms entering (6.4) with the simple roots of 5'L(n, M). 
Note that, in principle, the final set of dominant walls can only be decided 
when one starts from the complete list of all the dynamically relevant walls. 
In all the models we examine the set of dominant symmetry walls (6.4) will, 
however, be part of the final minimal set of dominant walls defining the 
complete billiard table. 
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If only the symmetry walls were present in the Hamiltonian, one would 
easily deduce from the above the following picture of the dynamics: When 
the trajectory hits a wall [3"''^^ = (3"' from within the interior of the billiard 
Ija+i — > it undergoes a reflection which reorders v""^^ and v"" from 
the incident state v"''^^ < v°- (in which — /3" decreases towards zero) 
into an outgoing state v"'^^ > v"- (in which /?"+^ — /3" increases away from 
zero). Each collision reorders a pair of velocity components so that, after a 
finite number of collisions, they will get reordered in a stable configuration 
where < v'^ ■ ■ ■ < v'^ and < (3"^ ■ ■ ■ < (3'^. The motion would then 
continue freely, i.e. without collisions, after the last reordering reflection. 
The same conclusion was already reached in Section 4 by a direct calculation 
of the dynamics of the Iwasawa scale factors of non-diagonal Kasner metrics 
(which is, indeed, where only symmetry walls are present). 

If we still consider for a moment the simple case of homogeneous but 
non-diagonal metrics and recall that the metric can be diagonalized at all 
times by a time-independent coordinate transformation — > x'* = i 
might appear that the symmetry walls, which are related to the off-diagonal 
components, are only a gauge artifact with no true physical content. This 
conclusion, however, would be incorrect. First, the transformation needed to 
diagonalize the metric may not be a globally well-defined coordinate transfor- 
mation if the spatial sections have non-trivial topology, e.g., are tori, since it 
would confiict in general with periodicity conditions. Second, even if the spa- 
tial sections are homeomorphic to M*^, the transformation a:* x'^ = U , 
although a diffeomorphism, is not a proper gauge transformation in the sense 
that it is generated by a non- vanishing charge, and therefore two solutions 
that differ by such a transformation should be regarded as physically distinct 
(although related by a symmetry). Initial conditions, for which the metric 
is diagonal and hence the symmetry walls are absent, form a set of measure 
zero. 

Let us finally mention the possibility of alternative treatments of the 
dynamics of off-diagonal metric components. We have just shown that the 
Iwasawa decomposition of the spatial metric leads to a projected description 
of the ^^((i, ]R)/S'0((i)-geodesics as motions in the space of the scale factors 
with exponential ( "Toda-like" ) potentials. An alternative description can 
be based on the decomposition G — AR of the spatial metric, where 
R e SO{d) and A is diagonal [92, 70]. One then gets Calogero-like potentials 
oc sinh~^(/5" — P'^). In the BKL limit, these potentials can be replaced by 
sharp wall potentials but whether the system lies to the left or to the right 
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of the wall /3" — = depends on the initial conditions in this alternative 
description. 

6.2 Curvature (gravitational) walls 

Next we analyze the gravitational potential, which requires a computation 
of curvature. To that end, one must explicitly express the spatial curvature 
in terms of the scale factors and the off-diagonal variables Again, the 

calculation is most easily done in the Iwasawa frame (4.8), in which the 
metric assumes the form (4.6). We use the short-hand notation = e'^"^ 
for the (Iwasawa) scale factors. Let C\J^x) be the structure functions of the 
Iwasawa basis {^"}, viz. 

dr^-i^V^Ar (6.5) 

where d is the spatial exterior differential. The structure functions obviously 
depend only on the off-diagonal components but not on the scale factors. 
Using the Cartan formulas for the connection one-form a;"^, 

dO" + ^uj\Ae^ = Q (6.6) 

b 

d^ab = i^ab + i^ba (6-7) 

where ujab = lac^^^b, and 

lab = SabAl = exp(-2/3")5„5 (6.8) 
is the metric in the frame {0°"}, one finds 

^ I r<^ j_ r"i d \ nb 

^ \ C C d r 




In the last bracket above, the commas denote the frame derivatives da = 
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The Riemann tensor R'^defi Ricci tensor R^^ a-nd the scalar curvature 
R are obtained through 

n\ = duJ\ + Y,uJ\^u\ (6.10) 

c 

= \Y.R\efO'^Of (6.11) 

where Jl"^ is the curvature 2-form and 

Rab = ^ -R'^acb' ^ X/ ^42 (6-12) 



a 



Direct, but somewhat cumbersome, computations yield 



where Fa is some complicated function of its arguments whose explicit form 
will not be needed here. The only property of Fa that will be of importance is 
that it is a polynomial of degree two in the derivatives d(3 and of degree one 
in 9^/3. Thus, the exponential dependence on the /3's which determines the 
asymptotic behaviour in the BKL limit, occurs only through the A^-terms 
written explicitly in (6.13). 

In (6.13) one obviously has h ^ c because the structure functions are 
antisymmetric in the pair [6c] . In addition to this restriction, we can assume, 
without loss of generality, that a 7^ 6, c in the first sum on the right-hand side 
of (6.13). Indeed, the terms with either a = b or a = c can be absorbed into 
a redefinition of Fa- We can thus write the gravitational potential density 
(of weight 2) as 

Vg ^ -gR = jE'^"'"""^''nCV)^ - E e-'^^(^^Fa (6.14) 

a,b,c a 

where the prime on ^ indicates that the sum is to be performed only over 
unequal indices, i.e. a ^ b,b ^ c,c ^ a, and where the linear forms CKabdP) 
and HaiP) are given by 

aabc{/3) = 2/3" + ^ {ay^b, by^c, cy^a) (6.15) 

e^a,b,c 
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and 



(6.16) 



respectively. Note that aahc is symmetric under the exchange of b with c, but 
that the index a plays a special role. 

Comparing the result (6.14) to the general form (5.25) we see that there 
are, a priori, two types of gravitational walls: the ct-type and the //-type. 
The a-type walls clearly come with positive prefactors, proportional to the 
square of a structure function . The /i-type terms seem to pose a problem 
because they do not have a definite sign. It would therefore seem that, in 
the BKL limit, the gravitational potential would tend to 



However, the indefinite //-type terms can generically be neglected in the BKL 
limit. This is most simply seen by noting that the inequalities aahc{P) > 
imply /ia(/3) > because /Xq is a linear combination with positive coefficients 
of the ttafec's- Indeed, we can write Hc = {ctabc + ctbca)/'^- Therefore the a-walls 
dominate the /t-ones. 

In fact, one can establish a stronger result, namely /laiP) > within 
the entire future light cone of the /5's. For this purpose, we note first that 
each linear form fJ^aiP) is lightlike, i.e. G^'^{^a)n{^J^a)u = 0. Therefore, each 
hyperplane /*„(/?) = is tangent to the light cone along some null generator. 
This means that the future light cone is entirely on one side of the hyperplane 
/^a(/3) = (i.e., either /(a(/3) > for all points inside the future light cone or 
/ia(/3) < 0). Now, the point = 13"^ = ■ ■ ■ = j3'^ = 1 is inside the future light 
cone and makes all the /i^'s positive. Hence f^aiP) > inside the future light 
cone for each a and 9[— 2/ia(/3)] = 0, and we really have 



Note, however, that the /t-type walls may make their existence felt in the 
exceptional case when (3 is close to the lightlike direction defined by /(. This 
is the case of "small oscillations" considered by BKL in [6, 8], for which they 
verify (for d = 3) that the evolution is indeed controlled by the ojabc-terms 
even in that region. 



lim Vg = ^'e[-2a„;.c(/3)] + E ( ± ^[-'^M]) (6.17) 





a,b,c 



(6.18) 
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Prom these considerations we deduce the additional constraints 

aabciP) > (i^ > 3) (6.19) 

besides the symmetry inequahtics (6.4). The hypcrplanes OiabciP) — are 
caUed the "curvature" or "gravitational" walls. Like the symmetry walls, 
they are timelike since 

G>'''{aabc)A<^abc). = +2 (6.20) 

The restriction D > 3 is due to the fact that in = 3 spacetime dimensions, 
the gravitational walls aabdP) — are absent, simply because one cannot 
find three distinct spatial indices. In this case all gravitational walls are of 
subdominant type and thus, in the BKL limit, 

VG^^(±e[-2/.„(/?)])^0 (D = 3). (6.21) 

a 

This is, of course, in agreement with expectations, because gravity in three 
spacetime dimensions has no propagating degrees of freedom (gravitational 
waves) . 

The fact that the gravitational potential becomes, in the BKL limit, a 
positive sum of sharp wall potentials, is remarkable for several reasons. First, 
the final form of the potential is quite simple, even though the curvature is 
a rather complicated function of the metric and its derivatives. Secondly, 
the limiting expression of the potential is positive, even though there are 
subdominant terms in Vq with indefinite sign. Thirdly, it is ultralocal in the 
scale factors, i.e. involves only the scale factors but not their derivatives. It 
is this fact that accounts for the decoupling of the various spatial points. 

The coefficients of the dominant exponentials involve only the undiffer- 
entiated structure functions C"bc- Consequently, one can model the grav- 
itational potential in leading order by spatially homogeneous cosmologies, 
which have constant structure functions, and by considering homogeneity 
groups that are "sufficiently non-abelian" so that none of the coefficients of 
the relevant exponentials vanish (Bianchi types VIII and IX for d — 3, other 
homogeneity groups for d > 3 - see [31]). By contrast, inspection of (6.13) 
and (6.14) reveals that the subleading terms not only lack manifest positivity, 
but also do depend on spatial inhomogeneities via the spatial gradient of the 
structure functions, and the first and second spatial derivatives of the scale 
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factors. It is quite intriguing that the associated walls are lightlike, unlike 
the walls associated with the leading terms all of which are timelike. Terms 
involving lightlike walls will thus have to be taken into account in higher 
orders of the BKL expansion. A Kac-Moody theoretic interpretation of this 
fact was recently proposed in [28]. 

The computations in this subsection involve only the Cartan formulas. 
They remain valid if in (4.8) we replace the one-forms dx^ by some anholo- 
nomic frame /* = f'''j{x)dx^ . This modifies the Iwasawa frame {6^}, which 
has no intrinsic geometrical meaning. The structure functions C^bdx) of the 
new frame get extra contributions from the spatial derivatives acting on f\ 
In fact, for p — dx'' not all gravitational walls 0;^^ appear because from 
(4.7) we then have C'^bc = for the top component (since 6'^ = dx'^), and 
C^^^bc = b,c d. Hence, the corresponding gravitational walls are absent. 
To get all the gravitational walls, one therefore needs an anholonomic frame 

However, the dominant gravitational wall aid-id is always present, and 
this is the one relevant for the billiard, when gravitational walls are relevant 
at all. 

6.3 p-form walls 

While none of the wall forms considered so far involved the dilatons, the 
electric and magnetic ones do as we shall now show. To make the notation 
less cumbersome we will omit the super- (or sub-) script (p) on the p-form 
fields in this subsection. 

6.3.1 Electric walls 

The electric potential density can be written as 

K) = ~\ e-^""i-»^'(^)(£"^-"'')2 (6.22) 

ai,a2,---,ap 

where are the components of the electric field tt'^ in the basis {^"} 

gai-ap = U'''j,U''^j^ ■ ■ ■ jp-K^^-^" (6.23) 

(recall our summation conventions for spatial cordinate indices) and where 
^ai-ap{P) are the electric wall forms 

ea,..a,(/9)=/3'^^ + --- + /3'^''-^0 (6.24) 
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Here the indices a/s are all distinct because 8°'^"'°'^ is completely antisym- 
metric. The variables S"'^'""''' do not depend on the j3^. It is thus rather easy 
to take the BKL limit. The exponentials in (6.22) are multiplied by positive 
factors which generically are different from zero. Thus, in the BKL limit, 
V^p-^ becomes 

K)^ E 0[-2ea....a,(/5)]. (6.25) 

<Il<02<'"<<Ip 

The transformation from the variables (jV^j, P'a, tt-^^ '■''') to the vari- 

ables {M'^i,V\,Aai-ap,£"'^"'"'^) is a point canonical transformation whose 
explicit form is obtained from 

a P 1 P a\,...,ap 

(6.26) 

The new momenta V'a conjugate to M'^i differ from the old ones P*a by 
terms involving M and A since the components Aa^-ap of the p-forms in 
the basis {^"} depend on the AA's, 

A = Af^^ ■ ■ ■ Af^p A ■ ■ 

■^ai---ap •'V ai •'v ap-^]v]p- 

However, it is easy to see that these extra terms do not affect the symmetry 
walls in the BKL limit. 



6.3.2 Magnetic walls 

The magnetic potential is dealt with similarly. Expressing it in the {^"}- 
frame, one obtains 

Vr^"^ E e-^'"«-->+i('^)(-^ar-a.+i)' (6-27) 



(P) 



2(p+l)! 

a\,a2,—,ap-\.\ 



where J^ai-ap+i are the components of the magnetic field Fmi-- mp+i in the 
basis {9"}, 

Ta....ap^, = ■ ■ ■M^^^\^^,F^,...i^^, (6.28) 

The ma^...ap^^{(S) are the magnetic linear forms 

6^{ai,a2,---ap+i} 
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where again all a^'s are distinct. One sometimes rewrites ma^...(jp_^^(/3) as 



m. 



ap+2— Od 



, where {ap+2, • • • , a^} is the set complementary to {oi, 02, • • • dp+i}', 



e.g., 

mi2-d-p-i = + • • • + P'^-P-' + y0 = (6.30) 

Of course, the components of the exterior derivative ^ of .4 in the non- 

holonomic frame {6"'} involves the structure coefficients, i.e. ^ai-op+i — 
c^[ai^a2 - ap+i] + C^-terms where da = J^\di is the frame derivative. 

Again, the BKL limit is quite simple and yields (assuming generic mag- 
netic fields) 

Kr^^ E ©[-2&a....a,_,_,(/3)]. (6.31) 

ai<---<ad-p-i 

Just as the off-diagonal variables, the electric and magnetic fields freeze in the 
BKL limit since the Hamiltonian no longer depends on the p-form variables. 
These drop out because one can rescale the coefficient of any ©-function to 
be one (when it is not zero), thereby absorbing the dependence on the p-form 
variables. 

The scale factors are therefore constrained by the further "biUiard" con- 
ditions 

ea,..„,(/3)>0, ma,...a,_,_M>0- (6.32) 

The hyperplanes ea-^...ap{f3) = and rhaj^...aj_^_j^{(3) = are called "electric" 
and "magnetic" walls, respectively. Both walls are timelike because their 
gradients are spacelike, with squared norm 



d-1 V 2 

(For D — 11 supergravity, we have d = 10,p = 3 and \p — and thus the 
norm is equal to +2). This equality explicitly shows the invariance of the 
norms of the p-form walls under electric-magnetic duality. 



6.4 Walls due to dilatons or to a cosmological constant 

The fact that the leading walls originating from the centrifugal, gravitational 
and p-form or are all timelike, is an important ingredient of the overall BKL 
picture. As we saw above, however, there exist subleading contributions in 
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the gravitational potential whose associated wall forms are lightlike. Such is 
the case also for the dilaton contribution 



(6.34) 



which has the same form as the subleading gravitational walls, since the 
exponentials that control its asymptotic behaviour are easily seen to be 
exp[— 2/ia(/3)]. Consequently, at least to leading order, we can neglect 
in the BKL limit. 

The only example of a spacelike wall that we know of is the cosmological 
constant term (in its weight-2 form) 



When A is positive (de Sitter sign) the spacelike wall (6.35) is repulsive. 
Depending on the initial conditions - which set the scale - this wall either 
prevents the system from reaching the BKL small volume regime or does not 
prevent the collapse. In the first case, the spacelike wall acts as a "barrier" to 
the motion of the billiard ball, and the reflection against it forces the billiard 
ball to run "backwards in /3-time" in the direction of increasing spatial volume 
(this is analogous to the bounce in the (global) de Sitter solution, which is a 
hyperboloid — a sphere that first contracts and then expands). In the second 
case, the billiard ball is already "beyond the barrier", and the presence of 
the spacelike wall has only a subdominant effect on its motion. In the BKL 
picture it quickly becomes neghgible: the cosmological potential h.g tends to 
zero as g goes to zero. When A is negative (Anti-de Sitter sign) the spacelike 
wall (3.37) is attractive and tends to favour collapse. It, however, quickly 
becomes negligible in the BKL limit. 

6.5 Chern- Simons and Yang-Mills couplings 

The addition of Chern-Simons terms, Yang-Mills or Chapline-Manton cou- 
plings does not bring in new (asymptotically relevant) walls. The only change 
in the asymptotic dynamics is a modification of the constraints. 

For Yang-Mills couplings, the contribution to the energy density from the 
Yang-Mills field takes the same form as for a collection of abelian 1-forms, 
with the replacement of the momenta tt* by the Yang-Mills momenta tt^ 
(where a — 1, - ■ ■ ,N and N is the dimension of the internal Lie algebra) and 




(6.35) 



a 



46 



of the magnetic fields by the corresponding non-abehan field strengths. As 
their abelian counterparts, these do not involve the scale factors f3^. Because 
of this key property, the same analysis goes through. Each electric and 
magnetic 1-form wall is simply repeated a number of times equal to the 
dimension of the Lie algebra. Gauss law is, however, modified and reads 

V,ni = V,< + f\y,A^ = 0. (6.36) 

Here, Vj is the standard metric covariant derivative. Similarly, the momen- 
tum constraints are modified and involves the non-abelian field strengths. 

The discussion of Chapline-Manton couplings or Chern- Simons terms pro- 
ceeds in the same way. The energy- density of the p-forms has the same de- 
pendence on the scale factors as in the absence of couplings, i.e., provides 
the same exponentials. The only difference is that the wall coefficients are 
different functions of the p-form canonical variables; but this difference is 
again washed out in the sharp wall limit, where the coefficients can be re- 
placed by one (provided they are different from zero). The momentum and 
Gauss constraints are genuinely different and impose different conditions on 
the initial data. 

Nevertheless, the Chern-Simons terms may play a more significant role in 
peculiar contexts when only specialized field configurations are considered. 
This occurs for instance in [58], where it is shown that the D = 11 super- 
gravity Chern-Simons term for spatially homogeneous metrics and magnetic 
fields may constrain some electromagnetic walls to disappear "accidentally" . 
This changes the finite volume billiard to one of infinite volume. 

7 Cosmological billiards 

Let us summarize our findings. The dynamics in the vicinity of a space- 
like singularity is governed by the scale factors, while the other variables 
(off-diagonal metric components, p-form fields) tend to become mere "spec- 
tators" which get asymptotically frozen. This simple result is most easily 
derived in terms of the hyperbolic polar coordinates (p, 7), and in the gauge 
(2.5). In this picture, the essential dynamics is carried by the angular vari- 
ables 7 which move on a fixed biUiard table, with cushions defined by the 
dominant walls ^^'(7)- However, it is often geometrically more illuminat- 
ing to "unproject" this billiard motion in the full Minkowski space of the 
extended scale factors P'^. In that picture, the asymptotic evolution of the 
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scale factors at each spatial point reduces to a zigzag of null straight lines 
w.r.t. the metric G^,ydj3^dj3^ . The straight segments of this motion are in- 
terrupted by collisions against the sharp walls 



(7.1) 



defined by the symmetry, gravitational and p-form potentials, respectively. 
As we showed all these walls are timelike, i.e. they have spacelike gradients: 



Indeed, the gradients of the symmetry and gravitational wall forms have 
squared norm equal to +2, independently of the dimension d. By contrast, 
the norms of the electric and magnetic gradients, which are likewise positive, 
depend on the model. As we saw, there also exist subdominant walls, which 
can be neglected as they are located "behind" the dominant walls. 

In the /3-space picture, the free motion before a collision is described by 
a null straight line of the type of (3.21), with the constraint (3.22). The 
effect of a collision on a particular wall w^(/?) is easily obtained by solving 
(5.31), or (5.32), with only one term in the sum. This dynamics is exactly 
integrable: it suffices to decompose the motion of the /9-particle into two 
(linear) components: (i) the component parallel to the (timelike) wall hyper- 
plane, and {ii) the orthogonal component. One easily finds that the parallel 
motion is left unperturbed by the presence of the wall, while the orthogonal 
motion suffers a (one-dimensional) reflection, with a change of the sign of the 
outgoing orthogonal velocity with respect to the ingoing one. The net effect 
of the collision on a certain wall w{l3) then is to change the ingoing velocity 
vector entering the ingoing free motion (3.21) into an outgoing velocity 
vector v'^ given (in any linear frame) by the usual formula for a geometric 
reflection in the hyperplane w{[3) = 0: 



Here, all scalar products, and index raisings, are done with the /9-space metric 
G^^. Note that the collision law (7.3) leaves invariant the (Minkowski) length 
of the vector v^. Because the dominant walls are timelike, the geometric 
reflections that the velocities undergo during a coUision, are elements of the 
orthochronous Lorentz group. Each reflection preserves the norm and the 
time-orientation; hence, the velocity vector remains null and future-oriented. 



(7.2) 




(7.3) 
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From this perspective, wc can also better understand the relevance of 
walls which are not timelike. Lightlike walls (like some of the subleading 
gravitational walls) can never cause reflections because in order to hit them 
the billiard ball would have to move at superluminal speeds in violation of the 
Hamiltonian constraint. The effect of spacelike walls (like the cosmological 
constant wall) is again different: they are either irrelevant (if they are "be- 
hind the motion"), or otherwise they reverse the time-orientation inducing a 
motion towards increasing spatial volume ("bounce"). 

The hyperbolic billiard is obtained from the /3-space picture by a radial 
projection onto the unit hyperboloid of the piecewise straight motion in the 
polywedge defined by the walls. The straight motion thereby becomes a 
geodesic motion on hyperbolic space. The "cushions" of the hyperbolic bil- 
liard table are the intersections of the hyperplanes (7.1) with the unit hy- 
perboloid, such that the billiard motion is constrained to be in the region 
defined as the intersection of the half-spaces wa{P) > with the unit hyper- 
boloid. As we already emphasized, not all walls are relevant since some of the 
inequalities wa{(^) > are implied by others [26]. Only the dominant wall 
forms, in terms of which all the other wall forms can be expressed as linear 
combinations with non-negative coefficients, are relevant for determining the 
billiard. Usually, these are the minimal symmetry walls and some of the p- 
form walls. The billiard region, as a subset of hyperbolic space, is in general 
non-compact because the cushions meet at infinity (i.e. at a cusp); in terms 
of the original scale factor variables (3, this means that the corresponding 
hyperplanes intersect on the lightcone. It is important that, even when the 
billiard is non-compact, the hyperbolic region can have finite volume. 

Given the action (2.1) with definite spacetime dimension, menu of fields 
and dilaton couplings, one can determine the relevant wall forms and compute 
the billiards. For generic initial conditions, we have the following results, as 
to which of the models (2.1) exhibit oscillatory behaviour (finite volume 
billiard) or Kasner-like behaviour (infinite volume billiard) 

• Pure gravity billiards have finite volume for spacetime dimension D < 
10 and infinite volume for spacetime dimension D > 11 [33]. This can 
be understood in terms of the underlying Kac- Moody algebra [27] : as 
shown there, the system is chaotic precisely if the underlying indefinite 
Kac-Moody algebra is hyperbolic. 

• The billiard of gravity coupled to a dilaton always has infinite volume, 
hence exhibits Kasner-like behavior [5, 2, 29]. 
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• If gravity is coupled to p-forms (with p 7^ and p < D — 2) with- 
out a dilaton the corresponding bilhard has a finite vohime [25]. The 
most prominent example in this class is D = 11 supergravity, whereas 
vacuum gravity in 11 dimensions is Kasner-like. The 3-form is crucial 
for closing the bilhard. Similarly, the Einstein- Maxwell system in four 
(in fact any number of) dimensions has a finite-volume billiard (see 
[60, 75, 99] for a discussion of four-dimensional homogeneous models 
with Maxwell fields exhibiting oscillatory behaviour). 

• The volume of the mixed Einstein-dilaton-p-form system depends on 
the dilaton couplings. For a given spacetime dimension D and a given 
menu of p-forms there exists a subcritical domain V in the space of 
the dilaton couplings, i.e. an open neighbourhood of the origin Ap = 
such that: (i) when the dilaton couplings \p belong to V the general 
behaviour is Kasner-like, but {ii) when the Ap do not belong to T) 
the behaviour is oscillatory [24, 29]. For all the superstring models, 
the dilaton couplings do not belong to the subcritical domain and the 
billiard has finite volume. Note, however, that the superstring dilaton 
couplings are precisely "critical", i.e. on the borderline between the 
subcritical and the overcritical domain. 

As a note of caution let us point out that some indicators of chaos must 
be used with care in general relativity, because of reparametrization invari- 
ance, and in particular redefinitions of the time coordinate; see [18, 54] for a 
discussion of the original Bianchi IX model. 

We next discuss the link between the various time coordinates used in 
the analysis. When working in the gauge (2.5), the basic time coordinate 
is T. Let us see how the other time coordinates depend on T. First, we 
note that the dynamical variable A is asymptotically a linear function of T. 
Indeed, its conjugate momentum tt^ < is asymptotically constant, so that 
the integration of dX/dT = —■^7i\ yields 



From this behavior we infer the time dependence of the intermediate time 
coordinate r by integrating its defining relation dr — p^dT — —{2/TTx)p'^dX — 



A = —^TTxT + const. 



(7.4) 



Hence p is an exponential function of T: 




(7.5) 
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—{2/'Kx)pdp. This yields 



r = -{l/'Kx)p^ + const. (7.6) 

At this stage, the asymptotic hnks between T, r, A and p are similar to the 
results derived above for a Kasner solution and does not depend on whether 
the system is chaotic or non-chaotic. 

The situation is more subtle for the proper time t in the chaotic case 
(in the non-chaotic case where one settles in a Kasner regime after a finite 
number of collisions, nothing is changed, of course). The proper time is 
obtained by integrating 

dt = ^dr = ^p'dT = -{2/nx)V9Pdp. (7.7) 

The main term in the integrand is ^/g = exp(— po"), where a = T^iY- The 
quantity a oscillates chaotically and is difficult to control. As argued in [69] , 
however, one nevertheless gets (in the four-dimensional case) the usual Kas- 
ner relation ^/g ~ t (and thus r ~ Ini and T ~ In | lni|) up to subdominant 
corrections. We refer the interested reader to [69] for the details. 

Concerning the frequency of collisions, note that the billiard picture 
makes it clear that the typical time interval between two collisions is con- 
stant as a function of T. In other words, the number of collisions goes like 
A I Inr ~ In I lni|. 

The hyperbolic billiard description of the (3+l)-dimensional homoge- 
neous Bianchi IX system was first worked out by Chitre [17] and Misner 
[84]. It was subsequently generalized to inhomogeneous metrics in [71, 56]. 
The extension to higher dimensions with perfect fluid sources was considered 
in [73] , without symmetry walls. Exterior p-form sources were investigated 
in [57, 58] for special classes of metric and p-form configurations. As far as 
we know, however, the uniform approach used in the present paper based on 
a systematic use of the Iwasawa decomposition of the spatial metric is new. 

8 Kac-Moody theoretic formulation 

Although the biUiard description holds for all systems governed by the ac- 
tion (2.1), the billiard in general has no notable regularity property. In 

particular, the dihedral angles between the faces, which can depend on the 
(continuous) dilaton couplings, need not be integer submultiples of tt. In 
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some instances, however, the bihiard can be identified with the fundamental 
Weyl chamber of a symmetrizable Kac-Moody (or KM) algebra of indefi- 
nite type^, with Lorentzian signature metric [26, 27, 23]. Such billiards are 
called "Kac-Moody billiards" . More specifically, in [26] , superstring models 
were considered and the rank 10 KM algebras Eio and BEiq were shown to 
emerge, in line with earlier conjectures made in [62, 64]^°. This result was 
further extended to pure gravity in any number of spacetime dimensions, for 
which the relevant KM algebra is AEd, and it was understood that chaos 
(finite volume of the billiard) is equivalent to hyperbolicity of the underlying 
Kac-Moody algebra [27]. For pure gravity in D = 4 the relevant algebra 
is the hyperbolic algebra AE^ first investigated in [38]. Further examples 
of emergence of Lorentzian Kac-Moody algebras, based on the models of 
[13, 21], are given in [23]. 

The main feature of the gravitational billiards that can be associated 
with KM algebras is that there exists a group theoretical interpretation of 
the billiard motion: the asymptotic BKL dynamics is equivalent (in a sense to 
be made precise below), at each spatial point, to the asymptotic dynamics of 
a one-dimensional nonlinear o"-model based on a certain infinite dimensional 
coset space G/K, where the KM group G and its maximal compact subgroup 
K depend on the specific model. As we have seen, the walls that determine 
the billiards are the dominant walls. For KM biUiards, they correspond 
to the simple roots of the KM algebra. Some of the subdominant walls 
also have an algebraic interpretation in terms of higher-height positive roots. 
This enables one to go beyond the BKL limit and to see the beginnings of 
a possible identification of the dynamics of the scale factors and of all the 
remaining variables with that of a non-linear cr-model defined on the cosets 
of the Kac-Moody group divided its maximal compact subgroup [28]. 

The KM theoretic reformulation will not only enable us to give a unified 
group theoretical derivation of the different types of walls discussed in the 
preceding section, but also shows that the (3-space of logarithmic scale fac- 
tors, in which the billard motion takes place, can be identified with the Cartan 
subalgebra of the underlying indefinite Kac-Moody algebra. The various types 
of walls can thus be understood directly as arising from the large field limit 
of the corresponding ci-models. It is the presence of gravity, which comes 

^Throughout this chapter, we wih use the abbreviations "KM" for "Kac-Moody", and 
"CSA" for Cartan subalgebra. 

^°Note that the Weyl groups of the £'-family have been discussed in a similar vein in 
the context of {/-duality [79, 89, 3]. 
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with a metric in scale-factor space of Lorentzian signature, which forces us 
to consider infinite dimensional groups if we want to recover all the walls 
found in our previous analysis, and this is the main reason we need the the- 
ory of KM algebras. For finite dimensional Lie algebras we obtain only a 
subset of the walls: one of the cushions of the associated biUiard is missing, 
and one always ends up with a monotonic Kasner-type behavior in the limit 
t ^ O"*" (for instance, the non-diagonal Kasner solution discussed in chapter 4 
has only symmetry walls corresponding to the finite dimensional coset space 
GL{d,M.) / SO{d)). The absence of chaotic oscillations for models based on 
finite dimensional Lie groups is consistent with the classical integrability of 
these models. While they remain formally integrable for infinite dimensional 
KM groups, one can understand the chaotic behavior as resulting from the 
projection of a motion in an infinite dimensional space onto a finite dimen- 
sional subspace. 

Before proceeding we should like to emphasize that the equivalence be- 
tween the models discussed in the foregoing sections and the KM cr-models 
to be presented in this section has so far only been established, in the case 
of general models, for their asymptotic dynamics. A proposal relating the 
infinitely many off-diagonal degrees of freedom arising in the KM cr-model 
to spatial gradients of the metric and other fields, as well as possibly other 
degrees of freedom has recently been made in [28]. The first steps of this 
proposal have been explicitly checked for the relation between D — 11 su- 
pergravity and the Eiq cr-model. The relevance of non-linear cr-models for 
uncovering the symmetries of M-theory has also been discussed from a dif- 
ferent, spacetime-covariant point of view in [100, 95, 96], but there it is En 
rather than Eiq that has been proposed as a fundamental symmetry. 

8.1 Some basic facts about KM algebras 

We first summarize some basic results from the theory of KM algebras, refer- 
ring the reader to [66, 85] for comprehensive treatments. As explained there, 
every KM algebra q = q{A) is defined by means of an integer- valued Cartan 
matrix A and a set of generators and relations [66, 85]. We shall assume that 
the Cartan matrix is symmetrizable since this is the case encountered for cos- 
mological billiards. The Cartan matrix can then be written as {i,j = 1, . . . r. 
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with r denoting the rank of £|(^)) 

_ 2{ai\aj) 
{ailai} 

where {ctj} is a set of r simple roots, and where the angular brackets denote 
the invariant symmetric bilinear form of q{A) [66]. Here the bilinear form 
acts on the roots, which are hnear forms on the Cartan subalgebra (CSA) 
f) C q{A). The generators, which are also referred to as Chevalley-Serre 
generators, consist of triples {hi, Cj, /j} with i = 1, . . . , r, and for each i form 
an s/(2,M) subalgebra. The CSA () is then spanned by the elements hi, so 
that 

[h^,hJ]=0 (8.2) 

Furthermore, 

[ci, fj] = Sijhj (8.3) 

and 

[hi, ej] = Aijej , [hi, /y] = -Aj/j (8.4) 

so that the value of the linear form aj, corresponding to the raising operator 
Cj, on the element hi of the preferred basis {hi} oi I) is aj{hi) — Aij. More 
abstractly, and independently of the choice of any basis in the CSA, the roots 
appear as eigenvalues of the adjoint action of any element h of the CSA on the 
raising (e,) or lowering (/j) generators: [h,ei] = +ai{h)ei, [h,fi] = -ai{h)fi. 
Last but not least we have the so-called Serre relations 

ad (e,)^-^- (e,) = , ad (Z,)^"^- (/,•) = (8.5) 

Every KM algebra possesses the triangular decomposition 

q{A) = n" © (} ® n+ (8.6) 

where n"^ and n~, respectively, are spanned by the multiple commutators of 
the Ci and /j which do not vanish on account of the Serre relations and the 
Jacobi identity. To be completely precise, n+ is the quotient of the free Lie 
algebra generated by the e^'s by the ideal generated by the Serre relations 
{idem for n~ and fi). In more mundane terms, when the algebra is realized, in 
a suitable basis, by infinite dimensional matrices, n"*" and n" simply consist 
of the "nilpotent" matrices with nonzero entries only above or below the 
diagonal. Exponentiating them formally, one obtains infinite dimensional 
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(8.1) 



matrices again with nonzero entries above or below the diagonal, and I's on 
the diagonal - as exemplified in the finite dimensional case by the matrices 
M in chapter 4. 

One of the main results of the theory is that, for positive definite A, 
one just recovers from these relations Cartan's list of finite dimensional Lie 
algebras (see e.g. [53] for a clear introduction). For non positive-definite A, 
on the other hand, the associated KM algebras are infinite dimensional. If 
A has only one zero eigenvalue (with all other eigenvalues strictly positive) 
one obtains the so-called affine algebras, whose structure and properties are 
rather well understood [66, 47]. For indefinite A [i.e. at least one negative 
eigenvalue of A), on the other hand, very little is known, and it remains an 
outstanding problem to find a manageable representation for them [66, 85] 
(see also [44] for a physicist's introduction). In particular, there is not a 
single example of an indefinite KM algebra for which the root multiplicities, 
i.e. the number of Lie algebra elements associated with a given root, are 
known in closed form. The scarcity of results is even more acute for the 
"Kac-Moody groups" obtained by formal exponentiation of the associated 
Lie algebras. In spite of these caveats, we will proceed formally, making sure 
that our definitions reduce to the standard formulas in the truncation to 
finite dimensional Lie groups, and more generally remain well defined when 
we restrict the number of degrees of freedom to any finite subset. 

In the remainder we will thus assume A to be Lorentzian, i.e. non- 
degenerate and indefinite, with one negative eigenvalue. We shall see that 
this choice is physically motivated by the fact that the negative eigenvalue 
can be associated with the conformal factor which, as we saw, is the one 
degree of freedom making the reduced Einstein action unbounded from be- 
low. As a special, and important case, this class of Lorentzian KM algebras 
includes hyperbolic KM algebras whose Cartan matrices are such that the 
deletion of any node from the Dynkin diagram leaves either a finite or an 
affine subalgebra, or a disjoint union of them. 

The "maximal compact" subalgebra fi is defined as the invariant subalge- 
bra of q{A) under the standard Chevalley involution, i.e. 

e(x) = X for all xet (8.7) 

with 

9{hi) = -hi , e{ei) = -fi , e{fi) = -Si (8.8) 

More explicitly, it is the subalgebra generated by multiple commutators of 
(ej — fi). For finite dimensional q{A), the inner product induced on the 
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maximal compact subalgebra 6 is negative-definite, and the orthogonal com- 
plement to i has a positive definite inner product. This is not so, however, 
for indefinite A (see the footnote on page 438 of [65]). 

It will be convenient in the following to introduce the operation of trans- 
position acting on any Lie algebra element E as 

E'^ := -e{E) (8.9) 

In this notation the relations above become hj = hi,ej = fi, ff = e^, and the 
subalgebra i is generated by the "anti-symmetric" elements satisfying E"^ — 
—E. After exponentiation, the elements of the corresponding maximally 
compact subgroup K formally appear as "orthogonal matrices" obeying — 

k-\ 

Sometimes it is convenient to use a so-called Cartan-Weyl basis for q{A). 
Using Greek indices /i., t/, . . . to label the root components corresponding to 
an arbitrary basis if^ in the CSA, with the usual summation convention and 
a Lorentzian metric G^,^ for an indefinite Q, we have hi :— 0(fHn, where 
ccf are the "contravariant components", G^j.^a'^ = cij^, of the simple roots 
ai {i = 1, . . .r), which are linear forms on the CSA, with usual "covariant 
components" defined as ai/^ = ai{Hn). 

To an arbitrary root a corresponds a set of Lie-algebra generators E^^s, 
where s — 1, . . . , mult (a) labels the (in general) multiple Lie-algebra ele- 
ments associated with a. The root multiplicity mult (a) is always one for 
finite dimensional Lie algebras, and is also one for the positive norm roots 

= {a\a) > (i.e. the "real roots") of general KM algebras (including, of 
course, the simple roots), but generically exhibits exponential growth as a 
function of —a^ for indefinite A. In this notation, the remaining Chevalley- 
Serre generators are given by := Ea^ and fi :— E^^i- Then we have 

[H^,E^,,]^af,E^,s (8.10) 

and 

[Ea,s,Ea',t] = ^ u (8-11) 
u 

The elements of the Cartan-Weyl basis are normalized such that 

{Hij,\H„) — , {Ea,s\E/3^t) — ^st^a+/3,o (8-12) 

where we have assumed that the basis satisfies E'^ — E^a- Let us finally 
recall that the Weyl group of a KM algebra is the discrete group generated 
by reflections in the hyperplanes orthogonal to the simple roots. 
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8.2 Decomposition of AE^ into SL{3) representations 



As a special example, let us consider the hyperbolic KM algebra AE3, with 
Cartan matrix^^ 



and work out the first few terms of its decomposition into 5'L(3,M) repre- 
sentations. This decomposition refers to the adjoint action of the s/(3,M) 
subalgebra defined below on the complete KM algebra. 

From the mathematical point of view, the algebra (8.13) was first studied 
in [38] as it is the simplest hyperbolic KM algebra containing a non-trivial 
affine subalgebra (see also [68] and references therein for more recent work). 
Prom the physical perspective this algebra is of special interest for pure grav- 
ity in four dimensions both because its regular subalgebras and its Chevalley- 
Serre generators can be physically identified with known symmetry groups 
arising in dimensional reductions of general relativity [86], and because its 
Weyl chamber is related to the original BKL billiard [27] (see footnote 12). 
First of all, the SL{2, M) subgroup corresponding to the third diagonal entry 
of Aij can be identified with the Ehlers group [37] which acts on solutions 
of Einstein's equations with one Killing vector. The affine subgroup corre- 
sponding to the submatrix 



is the Geroch group acting on solutions of Einstein's equations with two com- 
muting Killing vectors [42] (axisymmetric stationary, or colliding plane wave 
solutions, respectively), see e.g. [62, 63, 12, 87] and references therein. Both 
the Ehlers group and the Geroch group act in part by nonlocal transforma- 
tions on the metric (or vierbein) components. 

A crucial role, in our analysis, is played by the SL{3, M) subgroup gener- 
ated by (ci, /i, hi) and (62, /2, h2), corresponding to the submatrix 




(8.13) 




(8.14) 




(8.15) 



11 



This algebra is also called [38], H3 [66], HA[^^ [68] or Af'^ [27]. 
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This group can be realized as acting on the spatial components of the metric 
(or vierbein), as an extension of the so-called Matzner-Misner SL{2, M) group 
for solutions with two commuting Killing vectors. 

The billiard picture for pure gravity in four dimensions can be readily 
understood in terms of the Weyl group of AEs [27]. For S'L(3,M), which 
has two simple roots, the Weyl group is the permutation group on three 
objects. The two hyperplanes orthogonal to the simple roots of SL{3, M.) can 
be identified with the symmetry walls encountered in section 6.1. The third 
simple root extending (8.15) to the full rank 3 algebra (8.13), which "closes 
off" the billiard, can then be identified with curvature wall orthogonal to 
0^123 [27]. Readers may indeed check that the scalar products between the 
two symmetry wall forms and the curvature wall form reproduce the above 
Cartan matrix"*^^. 

Let us now expand the nilpotent subalgebra in terms of representations 
of the A2 = s/(3, M) subalgebra of AE^ exhibited in (8.15). For this purpose, 
given any root a, we define its sl{3, M) level i to be the number of times the 
root eta appears in it, to wit 

a — mai + na2 + (-o.^ (8.16) 

Note that this notion of level is different from the affine level (= m) which 
counts the number of appearances of the over-extended root cti [38]. The 
hyperbolic algebra is thus decomposed into an infinite tower of irreducible 
representations of its s/(3,]R) subalgebra. Such a decomposition is simpler 
than one in terms of the affine subalgebra, whose representation theory is 
far more complicated (and, unlike that of s/(3, M), only incompletely under- 
stood). 

After inclusion of the third Cartan generator h^, the level i — sector 
is just a gl{3,M) subalgebra with generators K^j (where i,j — 1,2,3) and 
commutation relations 

[K'j,K^] = S^K'i - 5\K^^ (8.17) 

Note that in the original analysis of [6, 8, 17, 84], the symmetry walls are not included; 
the KM algebra that arises has the 3x3 Cartan matrix 




and its fundamental Weyl chamber (radially projected on the hyperbolic plane H^) is the 
ideal equilateral triangle having its three vertices at infinity. 
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These Lie algebra elements will be seen to generate the GL(3, M) group act- 
ing on the spatial components of the vierbein. The restriction of the ^£^3- 
invariant bilinear form to the level-0 sector is 

{K'j \K\)^ SjS^ - 5)5^ (8. 18) 

The identification with the Che valley- Serre generators is 

ei = K\ , h^K\ , h^^K\-K\ 

h3 = -K\-K\ + KS (8.19) 

manifestly showing how the over-extended CSA generator enlarges the 
original s/(3,M) generated by (ei,/i,/ii) and (e2,/2,/i2) to the Lie algebra 
gl{3,W). The CSA generators are related to the "central charge" generator 
c by 

c=h2 + h = -K\ (8.20) 

This is indeed the expected result: in the reduction of gravity to two dimen- 
sions (i.e. with two commuting Killing vectors), the central charge does not 
act on the internal degrees of freedom, but as a scaling on the conformal 
factor [63, 12, 87] (here realized as the 1-1 component of the vierbein). The 
affine level counting operator d is given by [38] 

d = hi + h2 + h^ = -K\ (8.21) 

and in contradistinction to c it does act on the internal volume in the reduc- 
tion to two dimensions (i.e. the "dilaton" p in the notation of [12, 87]). The 
operator d may be viewed as the Lq operator of a full Witt-Virasoro algebra 
enlarging the Geroch group via a semidirect product [65]. 

The irreducible representations of SL{3, M) are most conveniently charac- 
terized by their Dynkin labels. Let us recall that these labels are non-negative 
integers which characterize any highest-weight representation of any finite- 
dimensional Lie algebra. Let us consider a representation with maximal 
vector v\, of (highest) weight A, i.e. such that for any element h in the 
CSA, h{vA) — A{h)vA, where the Lie algebra acts on the representation vec- 
tor space (or "module"), and where A is a linear form on the CSA. The 
Dynkin labels of this representation are defined as Pi{A) = 2{ai\A)/ {ai\ai), 
where ai are simple roots of the Lie algebra under consideration. Note that 
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in the case of simply laced algebras such as s/(3,]R), and AE^, the simple 
roots have all a squared length equal to 2, so that the definition of the labels 
reduces to Pi{A) = (ail A) (while the Cartan matrix similarly simplifies to 
Aij = {ai\aj)). In the case of s/(3,M) we have two simple roots, ai and 0:2, 
and therefore two Dynkin labels: pi(A) = (q;i|A), p2(A) = (0:2 1 A). In terms 
of the Young tableau description of sl{3, R) representations, the first Dynkin 
label pi counts the number of columns having two boxes, while p2 counts 
the number of columns having only one box. For instance, (^1,^2) = (IjO) 
labels an antisymmetric two-index tensor, while (pi,P2) = (0,2) denotes a 
symmetric two- index tensor. Note also that the dimension of the s/(3, M) rep- 
resentation described by the labels (^1,^2) is {pi + 1)(P2 + l)(Pi +P2 + 2)/2. 

Let us determine the representations of s/(3, M) appearing at the first level 
{i = 1) of AE^ At this level, one sees that, under the adjoint action of 
s/(3, M), i.e. of (ei, /i, hi) and (e2, /2, /i2), the extra Che valley- Serre generator 
fs is a maximal vector. Indeed, 

ei(/3) = [ei,/3] = 

e2(/3) = [e2,/3] = (8.22) 

The weight of fs is A — —a^ because, by definition /^(/a) = [h, fs] — 
—oi^{h)f^. Hence the Dynkin labels of the representation built on the maxi- 
mal vector /s are pi = —(ailas) = —Ai^ with i = 1,2. Explicitly, this gives 
(j*!? P2) = (0,2). These labels also show up as the eigenvalues of the actions of 
the (specially normalized) Cartan generators {hi, /12) on the maximal vector 
/a, namely: 

hi{f,)^[hi,fs] = 

h2{f3) = [h2,f3] = 2/3 (8.23) 

As we said above, the representation (pi,p2) = (0,2) corresponds to a sym- 
metric (two-index) tensor. 

Thus, at the levels ±1 we have AE^ generators which can be represented 
as symmetric tensors E^^ — E^^ and Fij = Fji. One verifies that all algebra 
relations are satisfied with (a(jj) = {uij + aji)/2) 

[K'j,E''^] = 5^jE'^ + 5]E^' 

[K',,Fki] = ~5 lF,i - 5\Fk, 

^^This analysis is analogous to the one performed in [38] for the affine subalgebra A^^' . 
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{F,,\E'^) = 8f5'] (8.24) 
and the identifications 

63 = ^;'=^ , /3 = F33 (8.25) 

As one proceeds to higher levels, the classification of S'L(3,M) represen- 
tations becomes rapidly more complicated due to the exponential increase 
in the number of representations with level I. Generally, the representations 
that can occur at level £ + 1 must be contained in the product of the level-£ 
representations with the level-one representation (0,2). Working out these 
products is elementary, but cumbersome. Moreover, many of the represen- 
tations constructed in this way will drop out (i.e. not appear as elements of 
the AE^ Lie algebra). The complications are not yet visible at low levels; 
for instance, the level-two generator E"-^^^'^ = e"'^'^Ei^^, with labels (1,2), is 
straightforwardly obtained by commuting two level-one elements 

[E^\ E^^] = e^^^'Ej^^ + e-^^^^Ej^^ (8.26) 

A more economical way to sift the relevant representations is to work out 
the relation between Dynkin labels and the associated highest weights, using 
the fact that the highest weights of the adjoint representation are the roots. 
More precisely, the maximal vectors being (as exemplified above at level 1) of 
the "lowering type" , the corresponding highest weights are negative roots, say 
A = —a with a of the form (8.16) with non-negative integers £, m, n. Working 
out the Dynkin labels of A = —a by computing, the scalar products with the 
two simple roots of s/(3, M) then yields 

Pi=p — n — 2m , p2 = 0.-21 -\-m — 2n (8.27) 

As indicated, we shall henceforth use the notation (pi,j52) = {p-,q) for the 
Dynkin labels. This formula is restrictive because all the integers entering it 
must be non-negative. 

Inverting this relation we get 

^ = ¥-Ip-¥ 

n = li-lp-lq (8.28) 

with n > 2m > 0. A further restriction derives from the fact that the highest 
weight must be a root of AE^, viz. its square must be smaller or equal to 2: 

A^^l{p^ + q^+pq-e)<2 (8.29) 



61 



Consequently, the representations occurring at level ^ must belong to the 
list of all the solutions of (8.28) which are such that the labels m,n,p,q are 
non-negative integers and the highest weight A is a root, i.e. < 2. These 
simple diophantine equations/inequalities can be easily evaluated by hand 
up to rather high levels. 

However, the task is not finished. Although the above procedure sub- 
stantially reduces the number of possibilities, it does not tell us how often 
(including zero times!) a given representation appears (i.e. its outer multi- 
plicity). For this purpose we have to make use of more detailed information 
about AEs, namely the root multiplicities computed in [38, 66]. Matching 
the combined weight diagrams with the root multiplicities listed in table 
on page 215 of [66], we obtain the following representations in the decompo- 
sition of AE3 w.r.t. its s/(3, M) subalgebra up to level i < 5: 



£^1 





(0,2) 






(1,2) 






(2,2) 






(1,1) 






(3,2) 






(1,3) 






(2,1) 


(occurs twice) 




(0,2) 






(1,0) 




(p, = 


(4,2) 






(2,3) 


(occurs twice) 




(0,4) 


(occurs twice) 




(3,1) 


(occurs three times) 




(1,2) 


(occurs four times) 




(2,0) 


(occurs three times) 




(0,1) 


(occurs twice) 



(8.30) 



(we have not listed the representations compatible with (8.28) and (8.29), 
which drop out, i.e. have outer multiplicity zero). Going to yet higher 
levels will require knowledge of AE^ root multiplicities beyond those listed 
, 66, 68]. There is, however, an infinite set of admissible Dynkin labels 



m 
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that can be readily identified by searching for "affine" highest weights for 
which m = in (8.16). They are 

(p, ?) = (£- 1,2) ^ (m,n) = (0,^-1) (8.31) 

Because the associated highest weights all obey = 2 all these representa- 
tions (and their "transposed" representations) appear with outer multiplicity 
one independently of I. A second series of affine representations is 

(p,g) = (^,0) ^ (m,n) = (0,£) (8.32) 

Now we have A^ = 0, and the corresponding representations have outer 
multiplicity zero because the corresponding states are already contained as 
lower weight states in the rperesentations (8.31). 

The Lie algebra elements corresponding to (8.31) are thus given by the 
two conjugate infinite towers of s/(3, M) representations (for £ = 1, 2, 3, . . .) 

, r^-'^-^ju (8.33) 

The affine representations are distinguished because they contain the affine 
subalgebra A^p C AE-^. This embedding was already studied in [38], but in 
view of future applications and because the identification is a little subtle 
we here spell out some more details. The affine subalgebra is identified by 
requiring its elements to commute with the central charge c (8.20). ^Prom 
this requirement we infer that we must truncate the full algebra to the sub- 
algebra generated by the elements K^j for i,j G {2,3} and those generators 
in (8.33) with ii = . . . i^^i = 1 and j, k e {2, 3}. In physicists' notation, the 
affine subalgebra is spanned by the current algebra generators {T~,T^, T+} 
and the central charge c. The affine level m e Z is the eigenvalue of the 
affine level counting operator d: [rf, T„J = mT„i. We therefore conclude that 
the affine level zero sector^^ consists of the generators of the Ehlers SL[2, M) 
group, viz. 

To" - /3 {= F33) , Ti = hs , = es (= E^') (8.34) 
At affine level m = — 1 we have 

TI, = 62 (=^'3) 

Tl, = [e2,e3]^2E''^ 

r+ ^ [[62, 631,63] = 2[E'\ E''] = 2E,'' (8.35) 

^''Recall that the affine level m must not be confused with the level i used in our 
decomposition of AE3, cf. the definition (8.16). 
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Similarly, at affine level m — +1 we get 

= [[/2,/3],/3] =2[F23,F33] =2Fl23 

Ti = 1/2, fs] = F23 

T+ = h {.= K\) (8.36) 

and so on for higher affine levels. The precise identification of the affine 
subalgebra is of crucial importance for understanding the embedding of the 
Geroch group into the full hyperbolic algebra. 

Other cases of interest, where similar decompositions can be worked out 
include the indefinite Kac- Moody algebras £'10 and En, both of which have 
been conjectured to appear in D = 11 supergravity and M Theory, see [62, 
64, 28], and [100], respectively. The first six rungs in the the s/(10,]R) level 
decomposition of Eiq have been worked out in [28], and will be extended to 
higher levels, as well as to En, in [39]. 

8.3 Nonlinear cr-Models in one dimension 

Notwithstanding the fact that we know even less about the "groups" asso- 
ciated with indefinite KM algebras we will formulate nonlinear cr-models in 
one time dimension and thereby provide an effective and unified description 
of the asymptotic BKL dynamics for several physically important models. 
The basic object of interest is a one-parameter dependent KM group ele- 
ment V = V(t), which is assumed to be an element of the coset space GjK, 
where G is the group obtained by formal exponentiation of the KM algebra 
0, and K its maximal compact subgroup, which is again obtained by for- 
mal exponentiation of the associated maximal compact subalgebra t defined 
above. [In the "transpose" notation defined above, the group K is the group 
of "orthogonal elements": A;^ = k~^\ For finite dimensional 0(A) our defini- 
tions reduce to the usual ones, whereas for indefinite KM algebras they are 
formal constructs to begin with. [Formal constructs similar to our objects V, 
VV~^, etc. have been used in somewhat different settings in [19, 20, 100].] 
However, to ensure that our definitions are meaningful operationally, we will 
make sure at every step that any finite truncation of the model is well defined 
and can be worked out explicitly in a finite number of steps. 

In physical terms, V can be thought of as an extension of the vielbein 
of general relativity, with G and K as generalizations of the G'L((i, M) and 
local Lorentz symmetries of general relativity. For infinite dimensional G, 
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the object V thus becomes an "oo-bein". It is then natural to associate to 
this vielbein a "metric" , viz. 

M := V^V (8.37) 

which is invariant under the left action ( V — > kV) of the "Lorentz group" K 
(actually the truncation of M. to the relevant GL{n, R) subgroup of the KM 

algebras entering the models we study turn out to correspond to the matrix 
defined by the contravariant components, g^^ , of the spatial metric used in 
section 4 above). Exploiting this invariance, we can formally bring V into a 
"triangular gauge" 

V = A-f/^M= M^A^M (8.38) 

where the abelian part A belongs to the exponentiation of the CSA, and the 
nilpotent part M to the exponentiation of n"'" thus recovering the formulas 
which we already used in sections 4.1 and 4.2. This formal Iwasawa decom- 
position, which is an infinite-dimensional generalization of the one we used 
before, can be made fully explicit by decomposing A and H in terms of bases 
of f) and n+: 

mult (a) 

M{t) = exp(^ ^c.,sit)Ea,) (8.39) 

where A_)_ denotes the set of positive roots. The components parametriz- 
ing a generic element in the CSA f), will turn out to be in direct correspon- 
dence with the logarithmic scale factors (5^ — (/3", 0) introduced in section 3. 
We anticipated this correspondence by using the same notation. [As ex- 
plained above, the apparently "wrong sign" of the exponents of A is due to 
the fact that 7W will correspond to the inverse of the spatial metric] The 
main technical difference with the kind of Iwasawa decompositions used in 
the foregoing sections is that now the matrix V{t) is infinite dimensional for 
indefinite q{A). Observe that for finite dimensional matrices, there was no 
need to worry about root multiplicities, as these are always one. By contrast, 
there are now infinitely many i/'s, and consequently M contains an infinite 
tower of new degrees of freedom. Next we define 

mult(a) 
a6A+ s=l 
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with 

3a,s = i^a,s + Vj> + l/i/l> + • • (8.41) 

(we put quotation marks to avoid having to write out the indices) . To define 
a Lagrangian we proceed in the usual way. First we consider the quantity 

mult(a) 

VV-^ = /3'^//^+ ^ ^■^vHma,sEa,s (8.42) 

which has values in the Lie algebra q{A). Here we have set 

«(/?) = (8.43) 

for the value of the root a { = linear form) on the CSA element (3 — P^H^. 
Next we define 

p := i (^VV-i + (VV-^)^) 

^ mult(a) 
aeA+ s=l 

where we arranged the basis so that E'^^ = E_a^s- Then we define a KM- 
invariant a-model action as J dtC where the Lagrangian is defined by using 
the KM-invariant bilinear form (.|.) (cf. (8.12)) 

C = \n-'{P\P) 

mult(a) 

= n-^ (-O'^/?'^ + - ^ e^v{Mma,s3a,s) (8.44) 

aeA+ s=l 

Here the Lorentzian metric 0^^, is the restriction of the invariant bilinear 
form to the CSA. It can be identified with the metric in the space of the scale 
factors, which is why we adopt the same notation. We have introduced in this 
Lagrangian a "lapse function" n (not to be confused with the lapse function 
N introduced before), which ensures that our formalism is invariant under 
reparametrizations of the time variable. Finally, we can simply describe our 
Lagrangian (8.44) as that of a null geodesic in the coset space G/K. 
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Taking the algebra AE^ as an example, this Lagrangian contains the Kas- 
ner Lagrangian (3.10) (without dilaton) as a special truncation. More specif- 
ically, retaining only the level zero fields (corresponding to a GI/(3, M)/0(3) 
cr-model) 

V{t)\^^^ = exp{h\{t)K\) (8.45) 

and defining from h"'i, a vielbein by matrix exponentiation = (exp/i)"^, 
and a corresponding contravariant metric = e^cC^ci one checks that the 
bilinear form (8.18) reproduces (half) the Lagrangian (3.10). [This compu- 
tation works more generally for any GL(n, R), and was already used in [28] 
in the GL(10, R) decomposition of Eiq.] The level-0 sector is thus associated 
with the Einsteinian dynamics of a spatial dreibein depending only on time. 

Let us then consider the fields (f)ij associated with the level-one generators 
E^^ . This leads to a truncation of our KM- invariant cx-model to the levels 
£ = 0,1, i.e. 

V(t)[_^ ^ = eMh\{t)K\) eM(l>abE'"') (8.46) 

In the gauge n = 1, the Lagrangian now has the form C ~ {g^^gY+g^^g^^<P<P-i 
where g denotes the covariant metric gij. As the 0ij's enter only through 
their time derivatives, their conjugate momenta are constants of the 
motion. Eliminating the 0's in favour of the constant momenta H by a partial 
Legendre transformation yields the reduced Lagrangian (for the dynamics of 
gij{t)) of the form C ~ {g^^gY ~ g^^ g^^^^- In other words, the elimination 
of the 0's has generated a potential = Vtp{g) (in the gauge n = 1) 

n(^) « +9ij9ki^"''^'' (8-47) 

It is then easy to check that this potential can be identified with the lead- 
ing (weight-2) gravitational potential in (6.14) (corresponding to the gauge 

N = 1), namely V^"""^"'*^ = (-(//?)i'=^'^'°g = +\ggaig^^ g^^C^^C' -j^ under the 
identification of the structure constants C^jk with the momenta conjugate to 

Y^ij ^ e^^^'C^\i (8.48) 

Note that the trace jk drops out of this relation, and that II*-' is of weight 
one, hke e^^^. 

Consequently, when neglecting the subleading gravitational walls oc F^in 
(6.14) and in four spacetime dimensions, the BKL dynamics at each given 
spatial point Xq is equivalent to the £ = 0, ±1 truncation of the A^^s-invariant 
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dynamics defined by (8.44). The fields (pij{t) parametrizing the components 
of the AE3 coset element along the £ = 1 generators are canonically conjugate 
to the structure constants C^jk(t,Xo). 

Similarly, one would like to associate the generators (8.33) with higher 
order spatial gradients. However, the proper physical interpretation of these 
fields as well as of the other higher level components remains yet to be found. 
In the case of the relation between supergravity in D = 11 and the Eiq coset 
model one could pursue the correspondence between spacetime fields and 
coset coordinates up to the (7/(10,M) level i = 3 included (corresponding 
to height 29) [28]. The correspondence worked thanks to several "miracu- 
lous" agreements between the numerical coefficients appearing in both La- 
grangians. 

Varying (8.44) w.r.t. the lapse function n gives rise to the constraint that 
the coset Lagrangian vanish. Let us define the canonical momenta 

^/^ ^= 7^ = ^''G.J" (8.49) 

and the (non-canonical) momentum-like variables 

SjO 1 

^a,s -J—= 0^"' (MP))ja,s (8.50) 

The latter variables are related to the momenta canonically conjugate to the 
Voc.s, i-e. Va,s ■= SC/dua,s, by expressions of the form Ila,s = Pa,s + "i^P + 
uup + ■ ■ •", see next subsection. 

In terms of these variables the KM-invariant Hamiltonian corresponding 
to the above Lagrangian reads 

^ mult (a) 

7i(/3,7r,i/,p) = n(-G''%7r,+ J] J] exp ( - 2a(/3))n„,,n,,,) (8.51) 

The constraint of a vanishing Lagrangian becomes that of a vanishing Hamil- 
tonian: 

7i fa (8.52) 

The momentum-like variables n„ ^ can be thought of as infinite-dimensional 
generalizations of the asymptotically frozen combinations P^N'^j used in 
section 5.1, cf. Eq.(5.14). They do not Poisson-commute with one another 
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in general. Instead one has Poisson brackets of the form {110, ,^, s} = 
^trpM^-i^t- Note that TT^ is timehke or null because the potential is mani- 
festly non-negative; the motion in the CSA is timelike — in fact, in the limit 
to be studied below, a broken future-oriented lightlike line. 

Because the coefficients of the exponentials in (8.51) are non-negative we 
can now apply exactly the same reasoning as in chapter 5. The crucial point 
is that the Cartan variables (3 do not enter the commutation relations of the 
n as indicated above. Therefore, as above, the off-diagonal components Va,s 
and the momentum-like variables Il.a,s get frozen asymptotically, provided 
all spacelike walls are "behind the motion" so that they do not conflict with 
the BKL limit. As before, we can introduce hyperbolic polar coordinates for 
parametrizing the dynamics of the Cartan variables — frj^. When taking 
the limit p oo, the Hamiltonian (8.52) behaves as 

n ~ HA+(/9,7r) = ^G^^TT^TT, + /r^exp [ - 2a(/5)] (8.53) 

aeA+ 

with constants > 0. This looks like a Toda Hamiltonian, except that 
the sum extends over all roots, rather than only the simple ones, and that 
the underlying KM algebra is indefinite, and not just finite or affine as in 
standard Toda theory. However, the dominant potential walls in the limit 
p —>■ oo are given by those terms Kaexp[ — 2q;(/3)] for which a is a simple 
root. Indeed, by definition, the non-simple roots a — riiai + n20i2 + • • • 
give rise to potential terms oc [exp(— 2Q;i(/3))]"i[exp(— 2q;2(/3))]"'^ • • • which 
are subdominant w.r.t. the set of potentials corresponding to the simple 
roots. Therefore, the cr-model Hamiltonian is asymptotically equivalent to 
the truncation of (8.53) to the simple roots: 

1 

n ~ 7isimpie(/9, tt) = -G^^TT^TT, + ^ exp [ - 2a,(/5)] (8.54) 

1=1 

Such hyperbohc Toda models, restricted to the simple roots, were first intro- 
duced and studied in [43]. 

Furthermore, as explained in subsection 5.2, the potentials become sharp 
wall potentials in the BKL limit — > oo. Making again the "genericity 

^'"'Perhaps a useful analogy is to think of a mountainscapc (defined by Toda exponential 
potentials for all roots) ; when the rnountaintops rise into the sky, only the nearest moun- 
tains (corresponding to the simple roots) remain visible to the observer in the valley (the 
Weyl chamber). 
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assumption" Ki > for the dominant simple root contributions (which is 
subject to the same caveats as before), we finally obtain the asymptotic 
Hamiltonian 

1 

HooiP, tt) := lim n(P, tt) = -n^n^ + J] ©( - 2ai(/3)) (8.55) 

i=l 

where the sum is over the simple roots only, and the motion of the f3^ is 
confined to the fundamental Wcyl chamber ««(/?) > 0. Note that, in the 
present KM setup, all the walls enter on the same footing. There is nothing 
left of the distinctions that entered our foregoing BKL-type studies between 
different types of walls (symmetry walls, gravitational walls, electric walls,...). 
The only important characteristic of a wall is its height ht a = ni + ^2 + • • • 
for a root decomposed along simple roots as a = niai + 7220:2 + ■ ■ ■• 

The sum in (8.52) not only ranges over the real roots, which give rise to 
infinitely many timelike walls, but also over null and purely imaginary roots 
of the KM algebra giving rise to lightlike and spacelike walls, respectively. In 
view of our discussion in section 7, the significance of the latter for the KM 
billiard remains to be fully elucidated. Here we only remark that imaginary 
roots are of no relevance for the Weyl group of the KM algebra, which by defi- 
nition consists only of refiections against the timelike hyperplanes orthogonal 
to real roots. 

To conclude: in the limit where one goes to infinity in the Cartan di- 
rections, the dynamics of the Cartan degrees of freedom of the coset model 
become equivalent to a billiard motion within the Weyl chamber, subject to 
the zero-energy constraint 7i{P,Ti) = 0. Therefore, in the cases where the 
cosmological billiards that we discussed in the first part of this paper are 
of KM-type, the gravitational models are asymptotically equivalent (modulo 
the imposition of the additional momentum and Gauss constraints) to the 
product over the spatial points of independent (1 + 0) -dimensional KM coset 
models G/K. 

8.4 Integrability, chaos and consistent truncations 

In this final subsection, we show that the one- dimensional KM cr-models 
are formally integrable. Then we address the issue of why this formal in- 
tegrability is not incompatible with the occurrence of chaos in the billiard 
description. We also discuss various possible finite-dimensional truncations 
of our infinite-dimensional KM-invariant cr-model. 



70 



The equations of motion of the off-diagonal fields constitute an infinite 
system of non-linear ordinary differential equations of second order. As usual, 
they are equivalent to the conservation of the g(74)-valued Noether charge 

mult(a) 

= J'^H^ + J2 Yl (J'-,sEa,s + J-a,sE-a,s) (8-56) 
aGA_|. s=l 

in the gauge n = 1. It is easily checked from the first line of this equation 
that 

J = M-'M (8.57) 
This equation can be formally solved as 

M{t) - M(0) • exp(i J) (8.58) 

For finite dimensional matrices M., the system is thus completely integrable. 
This fact was already used in section 4 where we wrote down the exact so- 
lution (4.1) of a ^^(n, ]R)/0(n) model. More precisely, the GL{n,M.) analog 
of the (Lie- algebra valued) conserved charge J is tt'^j = \g^^gkj, whose con- 
servation is clear from Eq. (3.13) . [Note again that the KM gauge n — 1 
corresponds to TV = 1 so that the KM time variable t corresponds to the 
gravitational time scale r.] In the finite-dimensional case we could derive 
the general solution by diagonalizing the constant matrix n^j, which led to 
the Kasner solution, written in terms of the eigenvalues f " of 7r*j and of the 
diagonalizing matrix L as in (4.1). Then, thanks to finite dimensionality 
of the diagonalizing matrix, we showed there that, after a finite transition 
time (subsequently interpreted as linked to the finite number of collisions 
on the symmetry walls needed to reorder the eigenvalues v""), the resulting 
asymptotic motion assumed a very simple monotonic form. 

By contrast, for infinite dimensional matrices the existence of the formal 
solution (8.58) does not necessarily imply regular behavior for the CSA de- 
grees of freedom in the asymptotic limit where all other degrees of freedom get 
frozen. First of all, and in marked contrast to the finite dimensional case, it is 
not known for indefinite KM algebras whether a generic Lie-algebra element 
J can always be "diagonalized" , i.e. conjugated into the CSA by an element 
of the KM group G, see however [67] . Indeed, if that were the case, it would 
mean that a generic solution of the KM coset dynamics could be conjugated 
to a simple monotonic and diagonal Kasner-like solution — vH -\- (3q in the 
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CSA, a conclusion difficult to reconcile with the chaotic motion that obtains 
in the BKL limit if the KM algebra is hyperbolic. Secondly, even if J could 
be conjugated into the CSA (as might be the case for certain non-generic el- 
ements) there remains the possibility that the required diagonalizing matrix 
(an infinite dimensional analog of the matrix L in (4.1)) introduces an infinite 
number of effective collisions in the time development of the analog of (4.1), 
because the number of effective collisions was found to increase with the rank 
of the matrix. In the latter case, the chaotic motion for the CSA degrees of 
freedom would be the result of the projection onto a finite-dimensional space 
of a regular geodesic motion taking place in an infinite-dimensional phase 
space. In the former case, it might be the geodesic motion on the infinite- 
dimensional coset-space G/K which could be intrinsically chaotic. We will 
leave further investigation of these delicate mathematical questions to future 
work. 

In order to better understand why the existence of an infinite number of 
conserved quantities {J^, Ja,s, J-a,s} does not rule out chaos, let us consider 
possible consistent truncations of our cr-model. By a "consistent truncation" 
we here mean a sub-model whose solutions are solutions of the full model. 
In order to formalize this notion, let us introduce a gradation V of root 
space: the X>-degree of a root a — n\a\ -\- 7120.2 -|- • • • is defined as T>{q.) : = 
n{Dx +712^2 + • • ■, where (X>i, I>2, • • •) is a given set of non-negative integers. 
Examples are the s/(3, M)-level i, the affine level m, or the height 'D^a) = 
ht (a) = Ui + n2 + ■ ■ ■■ In terms of an expansion according to increasing 
gradation V, the full cr-model Lagrangian (8.44) has the structure (with all 
coefficients suppressed and in the gauge n = 1): 

+ exp (2ai(/5)) [u']'' + exp (2^2 (/9)) + u'u'f 
+ exp (2a3(/3)) [u^ + u^u^ + u^u^ + v^u^v^f + • • • (8.59) 

The notation here is somewhat schematic: the lower indices on the roots 
and the upper indices on the off-diagonal fields refer to the gradation T> (so 

a;i,a2,--- are not simple roots here). Furthermore, the degree zero term 
coincides with the free CSA kinetic term only if T); > for all i; for other 
gradations, the level zero sector is described by the a-model Lagrangian for 
the respective level-0 subalgebra (which may be non-abelian as was the case 
for the level ell used in section 4.2). The various terms within parentheses 
correspond to the gradation of the ja^s above: ji ~ i'^,j2 ~ + u^u^jj^ ~ 
z>^ -|- z/^z>^ -|- z/^z>^ -|- v^v^i/^, .... The momenta canonically conjugate to the 
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off-diagonal variables i/^,!/'^, - ■ • have a similar graded structure 

Pi ~ exp (2ai(/3)) ji + cxp (2^2 J2^^' + exp (2a3(/3)) j3[i/2 + t/V^] + • • • 

= Hi + Hsz/^ + n3(z/2 + z/V^) + • • • , 
P2 ~ exp (2a2(/5))j2 + exp (2a3(/5))j3t/^ H 

= n2 + n3 + 

P3 ~ exp(2a3(/3))j3 + --- = n3 + ---, (8.60) 

where, consistently with the above definition, we have introduced the short- 
hand notation n„ = exp (2q;„(/9)) j„. Formally, the triangular relations above 
can be inverted iteratively to get infinite series of the form 

Hi ~ Pl+P2l^^ +P3W^ + , 

n2 ~ P2+P3V^ +Pa[v^ + + 

Ha ~ P3+P4i^^H (8.61) 

This yields a formal expression for the Hamiltonian in terms of the canonical 
variables (/3, tt; i/",p„), where n again refers to the gradation, 

7^(/?,7^;I/^p„) - 7r2 + exp(-2ai(/3))n? + exp(-2a2(/3))n2 + 

exp(-2a3(/3))n2 + ... (8.62) 

It is now easy to see that one can obtain a consistent finite- dimensional 
truncation of the dynamics by requiring that all canonical momenta above 
a certain V-degree uq vanish . For instance, if we require p„ = for n > 3 
implying II^ = for all n > 3 we get a consistent dynamics following from 
the finite Hamiltonian 

n^^\P,TT;u\p,,u^P2) ~ 7r2 + exp(-2«i(/3))[pi+p2^^f + 

exp(-2a;2(/3))P2 (8-63) 

The finite-dimensional Lagrangian corresponding to this Hamiltonian reads 

£(2) ^ ^2 ^ (2ai(/3)) [z>^]2 + exp {2a2{P)) [z>' + i^^i^^f (8.64) 

and is obtained from the original Lagrangian (8.59) by setting to zero all jn 
for n > 3, where, for instance, js = z>^ -|- z/^z>^ + z/^z>^ -|- z/^z/^z>^, and so on. 
Note that the highest-degree variable does not explicitly appear in the 
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truncated Hamiltonian 'H^'^\ This implies that the corresponding highest- 
degree canonical momentum p2 is a constant of the motion. However, non- 
trivial dynamics is obtained for the remaining degrees of freedom j3, tt; iy^,Pi- 
Note also that a particular set of solutions of the r'-degree-2 dynamics above 
is obtained when the constant of motion p2 happens to vanish. This particular 
case of the P-degree-2 dynamics is simply the "D-degree-l truncation of the 
original dynamics, obtained by setting to zero all momenta above level 2, i.e. 
P2 = P3 — ■ ■ ■ — 0. The corresponding Hamiltonian is simply 

7iW(/3,7r;i/\pi) ~7r' + exp(-2Q;i(/3))p? (8.65) 

Again the highest- degree momentum, which is now pi, is a constant of the 
motion, so that Ti^^^ directly defines the reduced dynamics of the Cartan 

variables j3, n. 

For V = ht, the truncated dynamics (8.65) at degree V = 1 defines the 
finite hyperbolic Toda model (8.54) involving only the simple roots. We have 
seen in the previous subsection that this height- 1 Hamiltonian universally de- 
scribes the asymptotic dynamics of the full o"-model. [Its BKL-limit directly 
yields the universal Weyl-chamber billiard (8.55).] On the other hand, when 
the chosen gradation is the s/(3) level of AE^, we get the low-level trunca- 
tions of subsection 8.2 above: the level-0 truncation (p^ — . . . — 0) yields 
the Kasner dynamics, while the level-1 dynamics describes the leading effect 
of the gravitational walls. A similar hierarchy of truncations appears in the 
£"10 model of [28], where the analysis included levels i — 0, 1,2,3 w.r.t. to 
the s/(10,M) = Ag (see also [25]). 

Let us now clarify the meaning of the infinite set of conserved charges 
J±a,s within the context of the finite- dimensional truncations which we just 
discussed. To this aim, one expresses J in terms of the (non-canonical) 
Hamiltonian variables {j3, tt, u, H): 

mult (a) 

+ E ^ exp[-2a(/3)]H„,,A/'-i£;„,,AA 

aeA+ s=l 
mult(a) 

+ I] I] n^,sJ^-^E_^,,J\f (8.66) 

aeA+ s=l 

Here only the third term spreads over the full algebra g = n~ © f) © n"'", 
while the first and second are entirely contained in the parabolic subalge- 
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bras f) © n"*" and n+, respectively. From this fact, it is easy to see that any 
truncation where the IIq, vanish for aU roots strictly above some degree, i.e. 
for I^(q;) > no, implies the vanishing of the correspondingly low "negative- 
root charges" J-a,s- As for the non-zero highest- degree momenta Ilao,s with 
Vi^ao) = no, Eq. (8.66) simply yields Ilao,s = J-ao,s- We thus recover the 
result above concerning the constancy of the highest-degree momenta. We 
can then move up (always in "D-degree) by adding simple roots cij to a^. 
Eq. (8.66) then implies Ha^^aiyS = J-ao+a„s + cUa^J^aoyS, where the constant 
c is determined by the commutations relations. Continuing in this manner 
until one reaches the degree zero, one obtains for all the negative-degree 
(non- vanishing) momenta expressions which are linear in the charges and 
polynomial of ascending order in the off-diagonal fields u^^s. Note, however, 
that these relations do not suffice to eliminate non-trivial degrees of freedom 
from the truncated Hamiltonian,i.e. to replace it by an (on-shell) equivalent 
Hamiltonian depending on fewer variables. [The only degrees of freedom 
that one can straightforwardly ehminate are the trivial highest-degree ones, 
whose associated "position variables" do not enter the Hamiltonian.] One 
might think that one would get further relations, and eventually enough re- 
lations to eliminate the off-diagonal degrees of freedom, by considering the 
higher-degree components of Eq. (8.66). However, it does not work this way. 
The degree-zero projection, namely 

mult (a) 

n^H^+J2 Yl ^a,s{■^''E_^,s^f)\ ^c^H, (8.67) 

(where is a constant vector) gives an expression for the degree-zero mo- 
menta TT^ which is linear in the charges and polynomial in the z/q^s's. On the 
other hand, the positive-degree projection of Eq. (8.66) yields for the positive- 
degree components of the charges Ja,s expressions which (for any positive 
degree) will involve "hidden" off-diagonal fields i'a,s of degree 'D{a) > no. 
These new variables did not enter the truncated Hamiltonians (8.63), which 
shows explicitly that the infinite number of conserved charges J^^^ of the 
original cr-model do not provide sufficiently many autonomous constants of 
the motion for the truncated Hamiltonian to guarantee its "integrability" in 
the trivial sense of allowing one to reduce the number of degrees of freedom 
to zero. 

Let us also dispose of another paradox concerning the formal integrability 
of KM-related models. Even if we consider the hyperbolic Toda model defined 
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by the truncation of our general cr-model to height one, i.e. 7isimpic(/3, tt) given 
by Eq. (8.54) with (strictly) positive constants Ki, one finds that it has formal 
integrability features. Indeed, the model (8.54) admits a Lax pair. Setting 
Si = \J (ai\ai) /2, rescaling the Chevalley-Serre generators as e[ = SiCi, // = 
Sifi, h[ = sffi, and defining Wij = {ai\aj), we replace the variables (/3'^,7r^) 
by (r,p*) by r = \/^exp[ — ai{f3)] and T,jWijp> = G'^''{ai)n'K^. One can 
then check that the Lie- algebra valued expressions 

M = J2\P'h'^ + ^'(^^ + ' ^ = - (8-68) 

i i 

satisfy the Lax-pair evolution equation M = [K, M] . One would then expect 
to be able to write down (infinitely) many conserved quantities associated 
to the isospectrality of the evolution of M(t) . [In the finite-dimensional case 
all the traces of the matrix M{t) are constants of the motion, see e.g. [90] 
and references therein.] However, in the KM case this formal Lax-pair inte- 
grability does not yield any concrete constants of the motion because there 
is no useful analog of the matrix trace in the KM algebra. The integrability 
properties of hyperbolic Toda theories have also been discussed in [43] . 

Let us note also in passing that there is full compatibility between our 
general result of the asymptotic constancy of all the off-diagonal degrees 
of freedom z/q, 11^ ., and the existence of an infinite number of conserved 
quantities Ja,s- Morally speaking the two infinite sets of asymptotic constants 
of integration {limi/Q s}, {limna s} correspond to the doubly infinite set of 
charges {J±a,s} (in non-zero degree), and one can check the compatibility of 
this correspondence by using the results above concerning truncated models. 
As for the zero-level conservation law (8.67) it is also compatible with the 
asymptotic limit because, in the gauge we use here (analogous to the (3- 
space picture of the billiard) all the components of the zero-degree momenta 
TT^ tend to zero asymptotically (because of the redhifts accumulated on the 
receding walls). 

In summary, there is thus no contradiction between the formal integra- 
bility of our cr-model and the chaos that appears in its generic asymptotic 
dynamics. This is mainly due to the presence of infinitely many degrees 
of freedom, and the dynamics of the finite-dimensional truncations of the 
(7-model of strictly positive degree. 
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9 Conclusions 



In this paper, we have shown that theories involving gravity admit a remark- 
able asymptotic description in the vicinity of a spacelike singularity in terms 
of billiards in hyperbolic space. Depending on whether the actual billiard has 
finite or infinite volume, the dynamical evolution of the local scale factors is 
chaotic (oscillatory) or monotonic (Kasner-like). The billiard, and in partic- 
ular its volume, is a fundamental characteristic of the theory, in the sense 
that it is determined solely by the field content and the parameters in the 
Lagrangian, and not by the initial conditions (in the generic case; i.e., there 
may be initial conditions for which some walls are absent - and the billiard is 
changed -, but these are exceptional). Although we have not investigated the 
physical implications of this property for cosmological scenarios (in particu- 
lar, string-inspired cosmologies [40, 14, 76, 98, 41]), nor its quantum analog, 
we believe that this result is interesting in its own right because it uncovers 
an intrinsic feature of gravitational theories. As discussed in section 8, the 
regularity properties of some billiards appear to give a powerful handle on 
possible hidden Kac-Moody symmetries, which remain to be exploited for 
their full worth. 
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A Asymptotic freezing: a simple model 

We have seen in the text that the "non-diagonal" phase-space variables (Q, P) 
(i.e. all variables except the extended scale factors parametrizing the 
diagonal-Iwasawa-part of the metric, and the dilaton), get frozen to constant 
values in the BKL limit. We provide here a more detailed understanding 
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of this property by discussing a simpler model which captures the essential 
features of the Hamiltonian (5.25). 

Consider a system with two canonically conjugate pairs {q,p), {Q, P) and 
time-dependent Hamiltonian 

H{q,p, Q, P; T) = + ^(/(Q, P)fp^e-P'^ {k being any real number) 

(A.l) 

where p is p = exp(T), with T the time variable. We have simplified here the 
Hamiltonian (5.26) by eliminating the variable A, i.e. by replacing it by an 
explicit function of the coordinate time T, in the approximation of a constant 
tta, so that A is a linear function of T (we scaled things so that X — T). One 
can think of {q,p) as mimicking the scale factors, while {Q,P) mimicks the 
off-diagonal components or the p-form variables. In (A.l), there is only one 
potential wall for q (namely, the second term). We shall consider later the 
case with several walls. 

Let us start by remarking the important fact that the variable f{Q, P) 
has zero Poisson bracket with H, therefore it is a constant of the motion. 
To simplify the following argument, we assume that we perform a canonical 
transformation such that P"^^"^ — f{Q^P). Dropping henceforth the label 
"new" , we get P — Pq where Pq is a constant which we assume to be different 
from zero. The basic aim of this Appendix is then to show that the (new) 
conjugate variable Q will also tend to a constant as T — > -|-oo. 

In the limit of large times, the motion in is a free motion interrupted 
by a collision against the potential wall, 

g = Po|r-ro| + go (A.2) 

where Tq is the time of the collision, go the turning point, and po the constant 
momentum of q far from the wall. The location of the turning point is 
determined by energy conservation (using the fact that p{Tq) vanishes before 
changing sign): 

PoVoe-''°^° = pI (A.3) 
The time scale AT of the collision (during which one feels the influence of the 
exponential potential) is roughly of the order l/(poPo): the later the collision, 
the sharper the wall. Let us evaluate the change in Q in the collision. To 
that end, we need to integrate Q = Pp'^ exp^—pq) over the collision, which 
yields 

/oo 
^2-pfce-p(Po|T-To|+9o) (A.4) 
-oo 
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The integrand is maximum ai T — Tq. We can approximate the integral by 
the value at the maximum times the time scale of the collision. Using (A. 3), 
one gets 

Hence, the variable Q receives a kick during the coUision (which can be of 
order one at early times), but the later the collision (i.e. the larger Tq), the 
smaller the kick. 

Assume now that there is another wall with the same prefactor and the 
same time dependence, say at q = d, so that q bounces between these two 
walls, 

additional 2 ^ 

At each collision Q receives a kick of order 1/po = e~^°. Because the 
speed of q remains constant (in the large T limit), the collisions are equally 
spaced in T. Therefore, the time of the nth collision grows linearly with n: 
roughly Tq""* ~ nd/po The total change in Q is then obtained by summing 
all the individual changes, which yields 

(Ag)Tota;~ J^e-"'^/^" (A.6) 

n 

This sum converges. Therefore, after a while, one can neglect the further 
change in Q, i.e., assume = 0. The Hamiltonian describing the large time 
limit is obtained by taking the sharp wall limit in the above H, and reads 
therefore 

H=^p' + Q{-q) + Q{q-d). (A.7) 

The pair {Q, P) drops out because it is asymptotically frozen. Our analysis 
justifies taking the sharp wall limit directly in H for this system, which is 
the procedure we followed in the text to get the gravitational billiards. 

Note that, had we studied a more complicated model with different pref- 
actors for the different walls, the prefactors /a(Q, P) would no longer have 
been exactly conserved. However, their non- conservation would only have 
been driven by the "far-away" walls, so that their time-variation would have 
been exponentially small. 



79 



References 



[1] E. Anderson, "Strong-coupled relativity without relativity," arXiv:gr- 
qc/0205118. 

[2] L. Andersson and A.D. Rcndall, "Quiescent cosmological singularities," 
Commun. Math. Phys. 218, 479-511 (2001) [arXiv:gr-qc/0001047]. 

[3] T. Banks, W. Fischler and L. Motl, "Dualities versus singularities," 
JHEP 9901 (1999) 019 [arXiv:hep-th/9811194]. 

[4] V. A. Belinskii, JET? Letters 56, 422 (1992). 

[5] V.A. Belinskii and I.M. Khalatnikov, "Effect of scalar and vector fields 
on the nature of the cosmological singularity," Sov. Phys. JETP 36, 
591-597 (1973). 

[6] V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, "Oscillatory ap- 
proach to a singular point in the relativistic cosmology," Adv. Phys. 19, 
525 (1970). 

[7] V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, "Construction of 
a general cosmological solution of the Einstein equation with a time 
singularity", Sov. Phys. JETP 35, 838-841 (1972). 

[8] V.A. Belinskii, I.M. Khalatnikov and E.M. Lifshitz, "A general solution 
of the Einstein equations with a time singularity," Adv. Phys. 31, 639 
(1982). 

[9] B. K. Berger, "Numerical Approaches to Spacetime Singularities," 
arXiv:gr-qc/0201056. 

[10] B. K. Berger, D. Garfinkle, J. Isenberg, V. Moncrief and M. Weaver, 

"The singularity in generic gravitational collapse is spacclikc, local, and 
oscillatory," Mod. Phys. Lett. A 13, 1565 (1998) [arXiv:gr-qc/9805063]. 

[11] E. Bergshoeff, M. de Roo, B. de Wit and P. van Nieuwenhuizen, "Ten- 
Dimensional Maxwell-Einstein Supergravity, Its Currents, And The Is- 
sue Of Its Auxihary Fields," Nucl. Phys. B 195, 97 (1982). 

[12] P. Breitenlohner and D. Maison, "On the Geroch group", Ann. Inst. 
Henri Poincare 46, 215 (1986). 



80 



[13] P. Breitenlohner, D. Maison and G. W. Gibbons, "Four-Dimensional 
Black Holes From Kaluza-Klein Theories," Commun. Math. Phys. 120, 
295 (1988). 

[14] A. Buonanno, T. Damour and G. Veneziano, "Pre- big bang bubbles 
from the gravitational instability of generic string vacua," Nucl. Phys. 
B 543, 275 (1999) [arXiv:hep-th/9806230]. 

[15] G. F. Chapline and N. S. Manton, "Unification Of Yang-Mills Theory 
And Supergravity In Ten-Dimensions," Phys. Lett. B 120, 105 (1983). 

[16] D. F. Chernoff and J. D. Barrow, Phys. Rev. Lett. 50, 134 (1983). 

[17] D.M. Chitrc, Ph. D. thesis, University of Maryland, 1972. 

[18] N. J. Cornish and J. J. Levin, "The mixmastcr universe is chaotic," 
Phys. Rev. Lett. 78, 998 (1997) [arXiv:gr-qc/9605029]; "The mixmaster 
universe: A chaotic Farey tale," Phys. Rev. D 55, 7489 (1997) [arXiv:gr- 
qc/9612066]. 

[19] E. Cremmer, B. Julia, H. Lu and C. N. Pope, "Dualisation of dualities. 
I," Nucl. Phys. B 523 (1998) 73 [arXiv:hcp-tli/9710119]. 

[20] E. Cremmer, B. Julia, H. Lu and C. N. Pope, "Dualisation of dualities. 
II: Twisted self-duality of doubled fields and superdualities," Nucl. Phys. 
B 535 (1998) 242 [arXiv:hep-th/9806106]. 

[21] E. Cremmer, B. Julia, H. Lu and C. N. Pope, "Higher-dimensional origin 
of D = 3 coset symmetries," arXiv:hep-th/9909099. 

[22] E. Cremmer, B. Julia and J. Scherk, "Supergravity Theory In 11 Di- 
mensions," Phys. Lett. B 76, 409 (1978). 

[23] T. Damour, S. de Buyl, M. Henneaux and C. Schomblond, "Einstein 
billiards and overextensions of finite-dimensional simple Lie algebras," 
arXiv:hep-th/0206125. 

[24] T. Damour and M. Henneaux, "Chaos in superstring cosmology," Phys. 
Rev. Lett. 85, 920 (2000) [arXiv:hep-th/0003139]; [See also short version 
in Gen. Rel. Grav. 32, 2339 (2000).] 



81 



[25] T. Damour and M. Henneaux, "Oscillatory behaviour in homogeneous 
string cosmology models," Phys. Lett. B 488, 108 (2000) [arXiv:hep- 
th/0006171]. 

[26] T. Damour and M. Henneaux, "E(10), BE(IO) and arithmetical chaos 
in superstring cosmology," Phys. Rev. Lett. 86, 4749 (2001) [arXiv:hep- 
th/0012172]. 

[27] T. Damour, M. Henneaux, B. Julia and H. Nicolai, "Hyperbolic Kac- 
Moody algebras and chaos in Kaluza-Klein models," Phys. Lett. B 509, 
323 (2001) [arXiv:hep-th/0103094]. 

[28] T. Damour, M. Henneaux and H. Nicolai, "i?io and a small tension 
expansion of M theory" Phys. Rev. Lett. 89, 221601 (2002) [arXiv:hep- 
th/0207267]. 

[29] T. Damour, M. Henneaux, A. D. Rendall and M. Weaver, "Kasner- 
like behaviour for subcritical Einstein-matter systems," arXiv:gr- 
qc/0202069, to appear in Ann. Inst. H. Poincare. 

[30] G. Dautcourt, "On the ultrarclativistic limit of general relativity," Acta 
Phys. Polon. B 29, 1047 (1998) [arXiv:gr-qc/9801093]. 

[31] J. Demaret, Y. Dc Rop and M. Henneaux, "Chaos In Nondiagonal Spa- 
tially Homogeneous Cosmological Models In Space-Time Dimensions 
> 10," Phys. Lett. B 211, 37 (1988). 

[32] J. Demaret, J.L. Hanquin, M. Henneaux, P. Spindel and A. Taormina, 
"The Fate Of The Mixmaster Behavior In Vacuum Inhomogcneous 
Kaluza-Klein Cosmological Models," Phys. Lett. B 175, 129 (1986). 

[33] J. Demaret, M. Henneaux and P. Spindel, "Nonoscillatory Behavior In 
Vacuum Kaluza-Klein Cosmologies," Phys. Lett. 164B, 27 (1985). 

[34] B. S. DeWitt, "Quantum Theory Of Gravity. 1. The Canonical Theory," 
Phys. Rev. 160, 1113 (1967). 

[35] D. Eardley, E. Liang and R. Sachs, "Velocity-Dominated Singularities 
In Irrotational Dust Cosmologies," J. Math. Phys. 13, 99 (1972). 

[36] A. Eskin and C. McMullcn, "Mixing, counting and equidistribution in 
Lie groups", Duke Math. J. 71, 181-209 (1993). 



82 



[37] J. Ehlers, Dissertation Hamburg University (1957) 

[38] A.J. Feingold and LB. Frenkel, "A hyperbolic Kac-Moody algebra and 
the theory of Siegel modular form of genus 2", Math. Ann. 263, 87 
(1983) 

[39] T. Fischbacher and H. Nicolai, "Low level representations for Eiq and 
£■11" , in preparation 

[40] M. Gasperini and G. Veneziano, "Pre - big bang in string cosmology," 
Astropart. Phys. 1, 317 (1993) [arXiv:hep-th/9211021]. 

[41] M. Gasperini and G. Veneziano, "The pre-big bang scenario in string 
cosmology," arXiv:hep-th/0207130. 

[42] R. Geroch, J. Math. Phys. 12, 918 (1971); 13, 394 (1972). 

[43] R.W. Gebert, T. Inami and S. Mizoguchi, Int. J. Mod. Phys. All, 5479 
(1996). 

[44] R.W. Gebert and H. Nicolai, "£"10 for beginners" , in Strings and Symme- 
tries, Proceedings of the Giirsey Memorial Conference I, Lecture Notes 
in Physics (Springer Verlag, 1994). 

[45] G. W. Gibbons, "Aspects Of Supergravity Theories," Print-85-0061 
(Cambridge) Three lectures given at GIFT Seminar on Theoretical 
Physics, San Feliu de Guixols, Spain, Jun 4-11, 1984- 

[46] G. Gibbons, K. Hashimoto and P. Yi, "Tachyon condensates, CarroUian 
contraction of Lorentz group, and fundamental strings," JHEP 0209, 
061 (2002) [arXiv:hep-th/0209034]. 

[47] P. Goddard and D.I. Ohve, "Kac-Moody and Virasoro algebras in rela- 
tion to quantum physics". Int. J. Mod . Phys. Al, 303 (1986) 

[48] S. W. Hawking and R. Penrose, "The Singularities Of Gravitational 
Collapse And Cosmology," Proc. Roy Soc. Lond. A 314, 529 (1970). 

[49] S. Helgason, "Differential Geometry, Lie Groups, and Symmetric Sc- 
paces" , Graduate Studies in Mathematics vol. 34, American Mathemat- 
ical Society, Providence 2001 



83 



[50] M. Henneaux, "Geometry Of Zero Signature Space-Times," Print-79- 
0606 (Princeton), published in Bull. Soc. Math. Belg. 31, 47 (1979) (note 
the misprints in the published version, absent in the preprint version). 

[51] M. Henneaux, M. Pilati and C. Teitelboim, "Explicit Solution For The 
Zero Signature (Strong Coupling) Limit Of The Propagation Amplitude 
In Quantum Gravity," Phys. Lett. B 110, 123 (1982). 

[52] R. E. Howe and C. C. Moore, "Asymptotic properties of unitary repre- 
sentations," J. Functional Analysis 32, 72-96 (1979) 

[53] J. E. Humphreys, "Introduction to Lie Algebras and Representation 
Theory" , Graduate Texts in Mathematics 9, Springer Verlag, 1980 

[54] G. Imponente and G. Montani, "On the Covariance of the Mixmaster 
Chaoticity," Phys. Rev. D 63, 103501 (2001) [arXiv:astro-ph/0102067]. 

[55] C. J. Isham, Proc. Roy. Soc. Lond. A 351, 209 (1976). 

[56] V. D. Ivashchuk, A. A. Kirillov and V. N. Melnikov, "Stochastic Proper- 
ties Of Multidimensional Cosmological Models Near A Singular Point," 
JETP Lett. 60, 235 (1994) [Pisma Zh. Eksp. Teor. Fiz. 60, 225 (1994)]. 

[57] V.D. Ivashchuk and V.N. Melnikov, "Billiard Representation For Multi- 
dimensional Cosmology With Multicomponent Perfect Fluid Near The 
Singularity," Class. Quantum Grav. 12, 809 (1995). 

[58] V. D. Ivashchuk and V. N. Melnikov, "Billiard representation for multi- 
dimensional cosmology with intersecting p-branes near the singularity," 
J. Math. Phys. 41, 6341 (2000) [arXiv:hep-th/9904077]. 

[59] V. D. Ivashchuk and V. N. Melnikov, "Exact solutions in multidimen- 
sional gravity with antisymmetric forms". Class. Quantum Grav. 18, 
R87 (2001) 

[60] R.T. Jantzen, "Finite-dimensional Einstcin-Maxwell-scalar-field sys- 
tem," Phys. Rev. D 33, 2121-2135 (1986). 

[61] R. T. Jantzen, "Spatially homogeneous dynamics: A unified picture," 
arXiv:gr-qc/0102035. 



84 



[62] B. Julia, LPTENS 80/16, Invited paper presented at Nuffield Gravity 
Workshop, Cambridge, Eng., Jun 22 - Jul 12, 1980. 

[63] B. Julia, in the Proceedings of the Johns Hopkins Workshop on Current 
Problems in Partcile Physics "Unified Theories and Beyond" (Johns 
Hopkins University, Baltimore, 1984). 

[64] B. Juha, in Lectures in Applied Mathematics, AMS-SIAM, vol 21 (1985), 
p.355. 

[65] B. Juha and H. Nicolai, "Conformal internal symmetry of 2d cr-models 
coupled to gravity and a dilaton", Nucl. Phys. 482, 431 (1996) 

[66] V.G. Kac, Infinite Dimensional Lie Algebras, 3rd edn., Cambridge Uni- 
versity Press, 1990 

[67] D. Peterson and V.G. Kac, Proc. Nat. Acad.Sci. USA 80, 1778 (1983) 

[68] Seek- Jin Kang, "Root multiplicities of the hyperbolic Kac-Moody alge- 
bra HA^^^\ J. Algebra 160, 623 (1994) 

[69] I. M. Khalatnikov, E. M. Lifshitz, K. M. Khanin, L. N. Shchur and 
Ya. G. Sinai, "On The Stochasticity In Relativistic Cosmology," J. Stat. 
Phys. 38, 97 (1985). 

[70] A. M. Khvedelidze and D. M. Mladenov, "Bianchi I cosmology and 
Euler-Calogero-Sutherland model," arXiv:gr-qc/0208037. 

[71] A. A. Kirillov, Sov. Phys. JETP 76, 355 (1993). 

[72] A. A. Kirillov and A. A. Kochnev, JETP Letters 46, 436 (1987). 

[73] A. A. Kirillov and V. N. Melnikov, "Dynamics Of Inhomogeneities Of 
Metric In The Vicinity Of A Singularity In Multidimensional Cosmol- 
ogy," Phys. Rev. D 52, 723 (1995) [gr-qc/9408004]; 

[74] A. A. Kirillov and G. V. Serebryakov, "Origin of a classical space 
in quantum cosmologies," Grav. Cosmol. 7, 211 (2001) [arXiv:hep- 
th/0012245]. 

[75] V.G. LeBlanc, "Asymptotic states of magnetic Bianchi I cosmologies," 
Class. Quantum Grav. 14, 2281-2301 (1997); "Bianchi II magnetic cos- 
mologies," Class. Quantum Grav. 15, 1607-1626 (1998). 



85 



[76] J. E. Lidsey, D. Wands and E. J. Copeland, "Superstring cosmology," 
Phys. Kept. 337, 343 (2000) [arXiv:hep-th/9909061]. 

[77] E. M. Lifshitz, I. M. Lifshitz and I. M. Khalatnikov, Sov. Phys. JETP 
32, 173 (1971). 

[78] U. Lindstrom and H. G. Svendsen, "A pedestrian approach to high en- 
ergy hmits of branes and other gravitational systems," Int. J. Mod. Phys. 
A 16, 1347 (2001) [arXiv:hep-th/0007101]. 

[79] H. Lu, C. N. Pope and K. S. Stelle, "Weyl Group Invariance and p-brane 
Multiplets," Nucl. Phys. B 476, 89 (1996) [arXiv:hep-th/9602140]. 

[80] J. M. Maldacena and C. Nunez, "Supergravity description of field the- 
ories on curved manifolds and a no go theorem," Int. J. Mod. Phys. A 
16, 822 (2001) [arXiv:hep-th/0007018]. 

[81] G. A. Margulis, "Applications of ergodic theory to the investigation of 
manifolds of negative curvature," Funct. Anal. Appl. 4, 335 (1969) 

[82] C.W. Misner, "Mixmaster universe," Phys. Rev. Lett. 22, 1071-1074 
(1969). 

[83] C. W. Misner, "Quantum Cosmology 1," Phys. Rev. 186, 1319 (1969); 
"Minisuperspacc," In *J R Klauder, Magic Without Magic*, San Fran- 
cisco 1972, 441-473. 

[84] C.W. Misner, in: D. Hobill et al. (Eds), Deterministic chaos in general 
relativity. Plenum, 1994, pp. 317-328 [gr-qc/9405068]. 

[85] R.V. Moody and A. Pianzola, Lie Algebras with Triangular Decompo- 
sition, Wiley, New York, 1995 

[86] H. Nicolai, "A Hyperbohc Lie Algebra Prom Supergravity," Phys. Lett. 
B 276, 333 (1992). 

[87] H. Nicolai, in "Recent Aspects of Quantum Fields" , Proceedings Schlad- 
ming 1991, Lecture Notes in Physics (Springer Verlag, 1991) 

[88] B. Piohne and A. Waldron, "Quantum cosmology and conformal invari- 
ance", [arXiv:hep-th/0209044] 



86 



[89] N. A. Obers, B. Pioline and E. Rabinovici, "M-theory and U-duality 
on T**d with gauge backgrounds," Nucl. Phys. B 525, 163 (1998) 
[arXiv:hep-th/9712084]. 

[90] D.I. Olive and N. Turok, "The Toda lattice field theory hierarchy and 
zero curvature conditions in Kac- Moody algebras", Nucl. Phys. B265, 
469 (1986) 

[91] H. Ringstrom, "The Bianchi IX attractor," Ann. H. Poincare 2, 405-500 

(2001) [arXiv:gr-qc/0006035]. 

[92] M. P. Ryan, "The oscillatory regime near the singularity in Bianchi-type 
IX universes," Ann. Phys. (N.Y.) 70, 301 (1972). 

[93] M. P. Ryan, "Hamiltonian cosmology," Spring er-Verlag, Heidelberg 
(1972). 

[94] M. P. Ryan and L. C. Shepley, "Homogeneous Relativistic Cosmolo- 
gies," Princeton, USA: Univ. Pr. ( 1975) 320 P. ( Princeton Series In 
Physics). 

[95] 1. Schnakenburg and P. C. West, "Kac-Moody symmetries of IIB super- 
gravity," Phys. Lett. B 517, 421 (2001) [arXiv:hep-th/0107181]. 

[96] I. Schnakenburg and P. C. West, "Massive IIA supergravity as a 
non-hnear reahsation," Phys. Lett. B 540, 137 (2002) [arXiv:hep- 
th/0204207]. 

[97] C. Teitelboim, "The Hamiltonian Structure Of Space-Time," PRINT- 
78-0682 (Princeton), in: General Relativity and Gravitation, vol 1, A. 
Held ed.. Plenum Press, 1980. 

[98] D. Wands, "String-inspired cosmology," Class. Quant. Grav. 19, 3403 

(2002) [arXiv:hep-th/0203107]. 

[99] M. Weaver, "Dynamics of magnetic Bianchi VIq cosmologies," Class. 
Quantum Grav. 17, 421-434 (2000) [arXiv:gr-qc/9909043]. 

[100] P. C. West, "^11 and M Theory", Class. Quant. Grav. 18, 4443-4460 
(2001), [arXiv:hep-th/0104081] 



87 



[101] R. Zimmer, Ergodic Theory and Semisimple Groups, (Birkhauser, 
Boston, 1984) 



88 



